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ABSTRACT 



The construction of an apparatus for the condensation 
performance testing of a horizontal bundle of four tubes with 
various refrigerants was completed. The apparatus was 
instrumented, and data reduction software was developed to 
provide bundle and single tube condensation data. 

Two tube bundles were tested, smooth copper tubes and low 
integral-fin copper-nickel tubes, with two refrigerants, R-114 
and R-113. An enhancement ratio of about 2.0 for the overall 
heat transfer coefficient was demonstrated for the finned 
tubes over the smooth tubes. 

Internal contamination, possibly due to a breakdown of the 
refrigerant molecules when subjected to high temperatures in 
the boiling chamber, inhibited further meaningful data 
collection. Recommendations for improvement of the test 
apparatus are made. 
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I. 



INTRODUCTION 



A. BACKGROUND 

With the advent of more complex shipboard weapons/combat 
systems and increasing fuel costs, the United States Navy has 
increasingly recognized a need for an energy-efficient, light- 
weight, and high-capacity air-conditioning system. Such an 
advanced system has been proposed and is currently in its test 
and evaluation phase for the DDG-51 program [Ref. 1] . 

Among the substantial differences in the advanced air- 
conditioning plant over extant systems is the use of titanium 
finned tubes in the condenser, in place of the smooth copper- 
nickel alloyed tubes previously used. The use of titanium in 
the air-conditioner's condenser presents two significant 
advantages. First, since the United States Navy uses ambient 
seawater as a heat sink, the corrosion resistance of titanium 
should measurably improve system integrity and reliability. 
Secondly, the advanced air-conditioning system expects to 
realize a weight savings per unit of 2286 kg (5040 lbs) and 
the major portion of this weight savings will come from the 
use of the lighter-weight titanium condenser tubes [Ref. 1]. 

Since the early 1970 's the United States Navy has used R- 
114 as its primary refrigerant. The advantages of R-114 over 
the more-widely-used commercial refrigerants in a naval 
application derive from its low toxicity, temperature 
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stability, stability when in contact with moisture, and its 
applicability to lower-pressure systems [Ref. 1] . R-114^ will 
be the refrigerant utilized in the advanced air-conditioning 
plant. 

As no substantial data base on the heat-transfer capacity 
exists for a system utilizing titanium condenser tube bundles 
configured as proposed in the advanced air conditioning plant 
and consequently, no way exists to measure this design's 
performance against competing designs, the need arose to 
develop a test apparatus to provide this data base in order to 
accurately predict future system performance. The design and 
subsequent construction of this test apparatus at the Naval 
Postgraduate School began in 1987 and is documented in the 
Theses of LCDR David S. Zebrowski [Ref. 2] and of LCDR Thomas 
J. Murphy [Ref. 3]. The test apparatus was also intended to 
serve as a test platform for advanced boiling surfaces 
proposed for use in the advanced air-conditioning plant. The 
incorporation of the advanced boiling surface test platform 
and a condenser tube test platform into a single test 
apparatus allows the widest latitude in examining the effects 
of various heat inputs, coolant flow rates, and various levels 
of refrigerant contaminants on individual components and 
overall system performance. Zebrowski and Murphy constructed 



^New refrigerants may have to be used in the future due to 
the "ozone problem." 
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the basic apparatus, and Murphy operated the system for 
preliminary evaporator measurements. 



B. OBJECTIVES 

The major objectives of this thesis were: 

1. Refine the design and complete the fabrication of the 
apparatus for the testing of condenser tube bundles 
utilizing R-114 as the working fluid. 

2. Develop and instrument a data-acquisition system with 
the associated software to provide a data base on the 
heat-transfer capacity of the condenser section of the 
test apparatus. 

3. Validate the test apparatus by comparing condenser 
performance against extant data bases derived from 
conventional condenser tube bundles. 



3 



II . 



LITERATURE SURVEY 



A. GENERAL OBSERVATIONS 

Condensation is the phase transformation process in which 
a vapor is transformed to liquid by removal of latent heat. 
The promotion of condensation in heat exchangers is used 
extensively in applications for propulsion engineering and 
air-conditioning/ refrigeration cycles. Due to its importance 
in the aforementioned fields, considerable research has been 
directed at the factors influencing the process and ways in 
which the process can be enhanced. Such factors as various 
modes of condensation, surface orientation to vapor flow, the 
shear forces exerted on condensate film, various external and 
internal surface enhancements, turbulent effects due to 
cascading condensate flow from another surface, and the 
effects of fluid properties on the process are detailed 
extensively in various reviews. This thesis deals exclusively 
with film condensation of refrigerants on small horizontal 
tube bundles in a quiescent vapor where tubes with external 
fins are compared to smooth tubes. 

B. SINGLE TUBE INVESTIGATIONS 

1 . Smooth Horizontal Tube Studies 

The first comprehensive condensation model was 
developed by Nusselt in 1916 [Ref. 4] based on the assumption 
that a quiescent vapor at saturation temperature coming into 
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contact with a wall surface below saturation temperature would 
condense and form a continuous film of condensate growing in 
thickness as the film flowed off the surface under the 
influence of gravity. The condensate film at the vapor-liquid 
interface would be at the vapor saturation temperature with a 
temperature gradient in the film down to the wall surface 
temperature. No radiation or convection would take place at 
the liquid-vapor interface as both liquid and vapor are at the 
same temperature, however the amount of vapor condensing 
corresponds to the quantity of heat flowing through the film 
by thermal conduction. In the case of a horizontal tube, he 
assumed that condensate flows in a laminar manner around the 
sides of the tube and off the tube in a continuous sheet. The 
heat transfer coefficient would be maximum at the top center 
of the tube decreasing gradually around the surface of the 
tube, as the inherent thermal resistance in the film grows 
with the thickness of the film, and eventually goes to zero at 
tube bottom. Nusselt's model is limited by disregarding 
surface tension forces that tend to hold up condensate at the 
tube bottom until overcome by gravity, resulting in the 
production of condensate droplets rather than a continuous 
sheet. Nusselt's expression for the average heat transfer 
coefficient for a single tube subjected to a constant heat 
flux is given by; 



a^= .655U^.pJ.g.Ahynf-D^'Q"]^/^ 



( 2 . 1 ) 
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Although constrained by the aforementioned limitations, 
Nusselt's model remains the conservative bench-mark against 
which all other models are compared. 

2 . Exterior Surface Enhanced Tubes 

Research into ways to improve condensation performance 
in condensers is motivated by the idea that any force that 
acts to thin the condensate film promotes an enhancement of 
heat transfer by minimizing the resistance to heat flow. 
General enhancement techniques include the fabrication of 
surfaces that promote dropwise condensation, wrapping exterior 
surfaces with wire, installation of porous drainage strips, 
and the fabrication of finned tubes with fins of various 
geometries and various spacings along the tube. This thesis 
deals with a low integral-fin tube and its heat transfer 
enhancement over a smooth tube. 

In considering a low integral-fin tube during 
condensation, there exists two distinguishable regions on the 
circumference of the tube; a flooded region and an un-flooded 
region. The flooded portion of the tube defines the 
condensate retention angle of the tube with respect to the 
tube circumference. The smaller the condensate retention 
angle of the tube, the larger the heat transfer capacity of 
the tube. The fin can be divided into three regions; the fin 
tips, the sides of the fins, and the fin root area. The 
majority of heat transfer takes place at the fin tips. 
Surface tension forces pull the condensate from the tips down 
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the fin sides into the flooded root area, where gravity drains 



the condensate to the bottom of the tube. The amount of 
condensate retained along the circumference of the tube is 
dependent upon the ratio of surface tension forces to gravity 
forces . 



comprehensive model for the prediction of the heat transfer 
coefficient of finned tubes based on their experimental 
results for various test fluids (including R-22) and finned 
tubes of various metallic compositions, and various fin 
geometries. Their model assumes gravity-dominated flow and 
neglects surface tension effects completely. Subsequent 
experimental results from other sources indicate that the 
Beatty and Katz model overpredicts the heat transfer 
coefficient as surface tension increases or as fin density 
increases. Nevertheless, the Beatty and Katz model conformed 
to their reported overall heat transfer enhancement of up to 
2.3 for R-22 on finned tubes compared to smooth tube data. The 
Beatty and Katz model is given by: 



Beatty and Katz [Ref. 5], in 1948, developed a 



= .689- [X^-pJ-g-Ahyn^-AT^] 



1/4. rj^] 1/4 

■■D ■' 



( 2 . 2 ) 



where. 



D 



eg 



1 , 1/4 




(2.3) 
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and, 



L = TT • (D^ - D^)/4 -D 
o r ' o 

A three-region theoretical model, based on their experimental 
results with three rectangular finned tubes, was presented by 
Karkhu and Borovkov [Ref. 6], in 1971, including surface 
tension forces. Measured heat-transfer coefficients 
demonstrated a 50 to 100 percent increase in vapor-side heat 
transfer coefficients for steam and R-113 condensing on finned 
tubes compared to smooth tubes. Unfortunately, they did not 
report enhancements separately for the two fluids. In 1980, 
results reported by Carnavos [Ref. 7] condensing refrigerants 
on various finned tubes demonstrated an enhancement of up to 
400 percent in the heat transfer coefficient over smooth tube 
results. Work done by Sauer and Williams [Ref. 8], in 1982, 
on the condensing performance of finned tubes with oil- 
contaminated R-113 demonstrated a serious degradation of 
performance when the surface tension to density ratio was 
large. The conclusion formed was that the higher-surface- 
tension oil remained in the fin gaps rendering the finned 
surface ineffective. Results reported by Honda et al. [Ref. 
9], in 1983, condensing R-113 on finned tubes of various 
geometries showed improvements of 900 percent for the vapor- 
side heat transfer coefficients at a constant vapor-to-tube- 
wall temperature difference. The results reported by Kabov 
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[Ref. 10] in 1984 with refrigerants R-12 and R-21, indicated 
that the bulk of the latent heat was in fact removed on the 
lateral fin surfaces and that an optimum fin height and 
spacing was dependent upon the ratio of surface tension to 
gravity forces. The results published by Masuda and Rose 
[Ref. 11] in 1985 in experiments with R-113 condensing on low 
integral-fin tubes, confirmed that the overall heat transfer 
coefficient increases, in general, with decreasing fin 
spacing. Their results showed a 600 percent increase in 
enhancement over smooth tube performance. More recent results 
by Marto et al . [Ref. 12] with R-113, demonstrated a 700 
percent enhancement in performance over a smooth tube and gave 
an optimum fin spacing of 0.5 mm for that fluid. Work 
performed and published in the same time frame by Sukhatme et 
al . [Ref. 13] with R-11 on conventional integral-fin tubes and 
special pyramid-shaped fin tubes, reported enhancement ratios 
of 5 to 7 for the low integral fin tubes and 10.3 to 12.3 for 
the pyramid-shaped fin tubes. 

C. TUBE BUNDLE INVESTIGATIONS 
1 . Smooth Tube Bundles 

In smooth tube bundles, two conflicting factors play 
a role in determining bundle performance. First, the 
condensate flowing from the tubes above a given tube in a 
bundle tends to thicken the condensate film on that tube, 
hence increasing the resistance to heat transfer. This effect 
is known as the condensate inundation effect. Secondly, 
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droplets from other tubes striking the film surface on a tube 
with a velocity provided by gravity or vapor flow, can create 
ripples or waves in the condensate film imparting a turbulence 
within the condensate film that produces an enhancement of 
heat transfer performance. 

In 1949, Jakob [Ref. 14] elaborating on Nusselt's 
model, predicted that the vapor-side heat transfer coefficient 
for a tube in a bundle, compared to the first tube in the 
bundle, was a function of that tube's relative position in the 
bundle. This model was based upon the assumption that a 
continuous laminar sheet of condensate flowing off the tube 
directly above the tube considered and striking the top of 
this tube further thickened the condensate film on this tube. 
Jakob's model is given, by: 

(2.4) 

“l 

In 1958, Kern [Ref. 15] proposed a model based on the 
assumption that discrete droplets or jets of fluid from other 
tubes caused ripples in the condensate film diminishing the 
inundation effect. Kern's model is less conservative than 
Nusselt's and, in many cases, remains the applied industrial 
design standard. Kern's model is given by: 




N 



- 1/6 



(2.5) 
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Work by Chen [Ref. 16], published in 1961, proposed a model 
that considered the momentum gain of falling condensate as 
well as the condensation of vapor on sub-cooled condensate 
droplets or sheets. Chen's model which is essentially 
Nusselt's model times a factor that incorporates the phase 
change number is given by: 

^ [1 + 2-Ph- (N-1) ] (2.6) 

“l 

where the Phase Change number is given by: 



Ph 



CpAT 



Equation (2.6) is valid provided. 



Ph- (N-1) < 2.0 



A model proposed by Eissenberg [Ref. 17] in 1972, assumes that 
condensate does not always drain in a vertical direction but 
may be diverted sideways due to local vapor flow conditions. 
In this case, the condensate strikes subsequent tubes on their 
sides rather than on their tops, minimizing inundation effects 
on the condensate film in the top portion of the tube. 
Eissenberg 's model is given by: 




.60 + 



.42- (N 



- 1 / 4 ) 



(2.7) 
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Enhanced Tube Bundles 



2 . 

Katz and Geist [Ref. 18] in 1948, studying six fin 
tubes in a vertical row, using R-12 among other working 
fluids, found that the effect of condensate inundation was 
over-predicted by Jakob's model and proposed that the exponent 
in Equation (2.4) be changed to 0.06 for finned tubes. 

A theoretical model was proposed by Honda et al. [Ref. 
19] in work published in 1987, for finned tube bundles that 
showed good agreement, within 5 to 7%, for data taken 
condensing acetone and R-12 on finned tubes. Their model 
relies upon solving a set of algebraic equations describing 
the vapor-to-coolant conjugate heat transfer problem. A 
compendium of the aforementioned models, both single tube and 
bundles, is provided by Marto [Ref. 20]. 

As of the date of this writing, the author is unaware 
of any comprehensive model, substantiated by experimental 
results, that accurately predicts heat transfer performance 
for enhanced tube bundles with varying fin geometries and 
varying pitch. 
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III. EXPERIMENTAL APPARATUS 



A. CONDENSER/BOILER TUBE BANK TEST PLATFORM 
1 . Overview 

The composite test platform, with associated support 
systems, is depicted schematically in Figure 3.1. The boiler/ 
condenser unit was fabricated from 6.35 mm thick, rolled 
stainless steel plates designed to withstand an absolute 
pressure of 3 08 kPa. The top cylinder which serves as the 
condenser for the apparatus is 1.30 m in length with an 
external diameter of 0.61 m. The effective condensing length 
is 1.22 m. The condensing chamber is capped on either end 
with circular stainless steel plates of 0.71 m diameter and 
connected to the chamber with a flange and stud assembly. 
System integrity is maintained with a 3.20 mm thick rubber 
gasket. The end caps support the nylon block condenser tube 
mounts, stainless steel backing plates, auxiliary condenser 
coils (coolant entrance side) , and mixing chambers (coolant 
exit side). The condenser chamber has five 12.7 mm thick 
Pyrex glass view ports backed with 12.7 mm thick Plexiglas. 
The view ports are 127 mm in diameter and are located axially 
along the main chamber, angled so as to provide top to bottom 
views of the test tubes at various locations along the 
effective condensing length. The condenser chamber and 
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Figture 3.1 System Schematic 




ancillary equipment are depicted graphically and in 
photographs in Figures 3 . 2 through 3.5. 

Internal to the condenser chamber, a stainless steel 
shroud is fitted over the effective condensing length of the 
instrumented test tubes (see Figure 3.6). The purpose of the 
shroud is to channel refrigerant vapor along the inside 
circumference of the shell, collecting vapor at the top, and 
then forcing vapor through the vertical axis of the 
instrumented tube bank into the shroud's well, where 
condensate is collected and the remaining vapor is condensed 
by an auxiliary condenser. The auxiliary condenser is 
composed of five helically wound copper tubes of 9.53 mm 
diameter suspended inside the shroud well by cantilevered rods 
welded to the entrance end cap. The stainless steel shroud 
was fabricated with a glass panel serving as one side of its 
stem to permit viewing of the instrumented tubes through the 
view ports. 

The condenser chamber is attached to the boiler 
chamber by a rolled stainless-steel cylinder nominally 280 mm 
in diameter and 2 03 mm in length, located mid-way along the 
condenser chamber's length allowing condensate to drain by 
gravity. A collar dam was fitted at the connecting flange, 
for the initial purpose of providing a condensate drain point 
to send condensate to the auxiliary storage unit. 
Subsequently, a faster, more-efficient method of transferring 
refrigerant was devised that allowed refrigerant vapor to be 



15 




6 



Figure 3.2 Photograph of Apparatus and Support Systems 
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Figure 3 . 4 Photograph of Coolant Entrance Endcap 




Figure 3.5 Photograph of Coolant Exit Endcap 
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l*4gure 3.6 



Cross-section of Condensing Chajnt>er Schematic 
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sent and condensed in the auxiliary storage unit. A 9.53 mm 
diameter copper tubing was provided for the vapor flow from 
the top of the condenser chamber to the storage tank. A 
similar tubing was provided between the storage tank and the 
boiling chamber to allow liquid flow via gravity. 

The boiling chamber is connected to the condenser unit 
via the interconnecting riser as previously discussed. The 
boiling chamber was also fabricated from rolled stainless- 
steel and is cylindrical with an outside diameter of 0.61 m 
and a length of 0.279 m. The boiling chamber is fitted with 
two Pyrex view ports, strengthened by Plexiglas plates, for 
observation during operation. The boiling unit is comprised 
of three groups of tubes (see Figure 3.3). In the simulation 
tube bundle, there are five boiling tubes each nominally rated 
at 1.5 kW and located at the bottom of the boiler unit. In 
the auxiliary tube bundle, there are four boiling tubes each 
nominally rated at 4 kW and located two on either side of the 
instrumented tube bundle. The instrumented tube bundle is 
located in the lower center of the boiling unit and is 
comprised of 35 tubes, with ten active tubes rated at one kW 
each, five instrumented tubes, and 20 dummy tubes. The power 
provided to these tubes is controlled by three variac 
controllers located on a console adjacent to the test platform 
and fed by a 208 volt power supply. The variac controllers 
are graduated in one percent increments of maximum power for 
each tube bundle. The exact location of each tube bundle is 
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apparent in Figures 3.2 through 3.5. It is a unique feature 
of this test platform that both boiling and condensation 
phenomena in a closed loop system can be evaluated 
simultaneously. Details on the design and construction of the 
basic test platform are available from Zebrowski [Ref. 2] and 
Murphy [Ref. 3]. Lightoff and securing procedures for the 
apparatus are listed in Appendix B. 

2 . Ancillary Systems 

Component equipment that support the R-114 tube bundle 
test platform and are located external to the apparatus 
include: (1) an R-114 storage and transfer system, (2) a 
condenser cooling and flow control system, (3) a coolant 
(62.4% by weight mixture of ethylene glycol and water) sump, 
and (4) an eight-ton refrigeration unit. 

The R-114 storage and transfer system, as previously 
described, consists of a stainless-steel cylindrical tank 
0.350 m in diameter and 0.91 m in length located on a rack 
above the coolant sump. Transfer of the refrigerant is 
accomplished by boiling in the main boiling chamber with vapor 
being sent to the storage tank via 9.53 mm diameter copper 
tubing located in the top center of the condenser chamber. 
The vapor is condensed in the storage cylinder by means of a 
helical copper coil, suspended the length of the cylinder on 
a cantilevered bar, that is kept cooled by the water-ethylene 
glycol mixture. Liquid refrigerant can be returned by gravity 
to the test platform from the bottom of the storage tank 
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through a 9.53 mm diameter copper tubing to the boiling 
chamber. Faster transfer of R-114 liquid was also possible if 
the system pressure was at or below the atmospheric pressure. 
Notice that the storage tank experiences an absolute pressure 
of about 210 kPa. 

Coolant and flow control to the test apparatus and 
ancillary equipment is accomplished by two different flow path 
systems. Both flow systems are driven by two 0.5 HP constant- 
speed pumps that take suction on the main coolant sump. The 
coolant for the test condenser tube bank passes from the pump 
discharge through 76 mm diameter PVC piping to a Plexiglas 
header. At the header, flow is split and proceeds through 
15.9 mm diameter Tygon flexible tubing to a bank of 
rotameters. Flow can be controlled by throttling a gate valve 
located at the entrance to each flow meter. Coolant flow 
leaves the exit of each flow meter and proceeds to the 
instrumented condenser tube bank through flexible Tygon 
tubing. At the exit of the main condenser chamber, coolant 
leaves each tube and flows through flexible Tygon tubing to 
individual mixing chambers (Figure 3.7). From the mixing 
chambers, coolant passes through flexible Tygon tubing to a 
central Plexiglas header suspended above the test platform, 
and the collective coolant is piped back to the main sump 
through 76 iron diameter PVC pipe. Auxiliary coolant flow 
leaves the pump discharge through PVC piping and proceeds to 
a two-way ball valve, where the flow is either sent to the 
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Figure 3.7 Mixing Chamber Schematic 
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R-114 storage tank and back to the main sump or is sent 
through PVC pipe to a large rotameter. From the auxiliary 
flowmeter, coolant enters a central header and is then split 
to flow to the five auxiliary condenser coils and back to the 
coolant sump. Flow through each auxiliary condenser coil is 
controlled by a valve located at the coil inlet. The coolant 
sump has a 1.81 cubic meter capacity and is constructed of 
12.7 mm thick sheets of Plexiglas. The coolant is 
approximately 62.4% ethylene glycol and 37.6% water (by 
weight) . 

The coolant is chilled by an externally-located, 8 ton 
refrigeration system that continually re-circulates sump 
coolant with a 0.75 HP pump. The refrigeration system is 
capable of maintaining a sump temperature between -21 C and 
ambient temperature. 

3 . Instrumentation 

The coolant temperature rise is measured by series- 
connected thermopiles having ten junctions on either end. 
These thermopiles were fabricated using type-T thermocouple 
wire. The ends of the thermopiles are inserted in stainless 
steel wells with copper-plugged tips located at the entrance 
to each individual condenser tube and at the exit of this 
tube's mixing chamber. Inlet coolant temperature is measured 
with type-T thermocouples also located in the entrance wells. 
Refrigerant liquid and vapor temperatures are measured with 
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type-T thermocouples inserted in stainless-steel wells at the 
locations indicated in Figures 3.2 through 3.5. 



System pressure is monitored through a calibrated 
pressure-vacuum gage valid over a range of 30 inches mercury 
(vacuum) to 30 pounds per square inch (gage pressure) . 

4 . Tube Bundle Data Acquisition and Reduction 

A Hewlett-Packard 9816A computer was used to control 
a Hewlett-Packard 3497A Automatic Data Acquisition System, 
which read the output of the thermopiles and thermocouples. 
Readings were made in millivolts and were converted to 
temperature readings in the data reduction program. The 
channels read by the data acquisition system are listed in 
Table 3.1. 

5 . Tubes Tested 

Two sets of tubes were tested. The first was a smooth 
copper tube (inside diameter 13.26 mm, outside diameter 15.88 
mm) and the second was a low integral-fin copper-nickel tube 
(inside diameter 10.16 mm, root diameter 14.00 mm, outside 
diameter 15.88). These two sets of tubes were tested, in 
bundles and individually, with R-114 and R-113. 
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TABLE 3.1 



CHANNEL ASSIGNMENTS ON DATA ACQUISITION SYSTEM 



Channel Numbers 
0-2 
3-4 

5 

6 

7 

8 

20 

21 

22 

23 



Measurement 

Vapor Saturation Temperature 
Liquid Temperature 



Inlet 


Temperature 


coolant 


1st 


Tube 


Inlet 


Temperature 


coolant 


2nd 


Tube 


Inlet 


Temperature 


coolant 


3rd 


Tube 


Inlet 


Temperature 


coolant 


4th 


Tube 



Thermopile 1st Tube 
Thermopile 2nd Tube 
Thermopile 3rd Tube 
Thermopile 4th Tube 
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IV. SYSTEM OPERATION AND DATA REDUCTION 



A. TEST PLAN AND MODIFICATION OF DESIGN FOR CONDENSER 

TEST APPARATUS 

1 . Ability to Change Working Fluids 

Although envisioned in the original apparatus design 
by Zebrowski [Ref. 2] and Murphy [Ref. 3], no physical system 
existed for the change-out and storage of working fluids. The 
system was designed as a general test platform to evaluate the 
performance of various refrigerants and steam during 
condensation on a variety of test tubes. 

The R-114 storage and transfer system as described in 
Chapter III, ancillary systems, was designed and built to 
allow for the storage of refrigerants that evaporate at 
atmospheric temperature and pressure. This additional 
capability not only conserves expensive refrigerants but 
facilitates speedy change of these working fluids during the 
evaluation of a particular tube with a variety of 
refrigerants . 

2 . Tube Alignment and Reduction of Tubes in Test Bundle 

Visual inspection, during preliminary experimental 

runs, revealed that bowing of the tube bundle was 
detrimentally affecting condensation patterns. As a result of 
misalignment of the endcaps during fabrication, bending 
moments of varying magnitudes and directions resulted in 
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condensate flow striking subsequent tubes in the vertically- 
oriented bundle at different positions around the 
circumference of the tubes. Upon disassembly of the end 
plates of the apparatus and comparison with the proposed 
drawings for the DDG-51 refrigeration condenser provided by 
the David Taylor Research Center, miscalculations in the tube 
pitch were revealed. The vertical pitch required for the 
project is 35.74 mm, centerline to centerline. The nylon 
bundle plates, as manufactured, had a pitch that was 
considerably less. In an effort to correct the tube mis- 
alignment problem, a low-powered laser was proposed to align 
the bundle tubes. A new nylon endplate was fabricated for the 
exit endcap, at the correct pitch with the tube bundle reduced 
to four vertical tubes, due to constraints in the endcap 
openings. An aluminum mount for the laser was fabricated by 
machine shop personnel that fit snugly into the tube 
penetrations of the nylon block. Several test projections 
were conducted on a plastic template fitted to the entrance 
endcap to minimize misalignment possibilities. The plastic 
template was scribed and cut with particular attention paid to 
correct pitch and used to manufacture the entrance nylon 
endplate. Subsequent visual inspection during experimental 
runs revealed that misalignments and bowing had been 
corrected. 
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3 . Vapor Superheat Problems 



During initial experimental runs, observed 
discrepancies in the measured temperature increases across the 
condenser tube lengths when compared to predicted temperature 
increases, prompted consideration of the possibility of vapor 
superheat occurring as vapor (at saturation temperature less 
than room temperature) flowed from the boiling chamber through 
the riser and up the sides of the condensing chamber. A two- 
pronged strategy was developed to minimize this possible 
effect. First, the apparatus shell was completely insulated 
with 12.7 mm thick foam insulation and all ancillary tygon 
tubing was insulated with double-wrapped, 3.18 mm thick foam 
insulation. Secondly, all experimental runs were performed at 
vapor saturation temperatures above the ambient temperature. 

4 . Contamination Problems 

Visual inspection of the condensing tubes during 
operation, correlated with observed disparities in measured 
temperature increases across condenser tubes, revealed a 
structured, crystalline surface or matrix formation on the 
tube surfaces that inhibited heat transfer. The first 
appearance of the inhibiting matrix occurred during runs with 
smooth copper tubes and with the coolant inlet temperature at 
approximately -20 C. The working fluid, R-114, was 
transferred to the storage unit except for an oily residue 
that remained in the bottom of the boiling chamber. This 
residue was not analyzed, but was assumed to be machining oil 
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that remained after apparatus fabrication. The system was 
flushed twice, once with acetone and once with R-113, and the 
boiling chamber was scrubbed out. The smooth copper tubes 
were re-installed after a complete cleaning of their outside 
surfaces, and the experimental runs were repeated under the 



same conditions. 


Visual inspection 


again revealed 


the 


presence of the 


matrix 


formation on 


the tubes , and 


the 


possibility of 


water 


contamination 


in the R-114 


was 



conjectured as a probable cause. The matrix was thought to be 
ice crystals solidified on the tube surface. The decision was 
made to continue experimental runs with a different 
refrigerant (R-113) and a different type of condensing tube 
(copper-nickel 1024 fpm low integral-fin tube) , with the aim 
of understanding the conditions that caused the phenomena to 
exist. First, a run was made using R-114 with the copper- 
nickel finned tube under exactly the same conditions as 
described in the smooth copper tube runs. Visual inspection, 
during this run, revealed an apparent thickening of the fins 
at the top of the tube when compared visually with the second 
tube (see Figure 4.1(a) fin normal appearance and Figure 
4.1(b) fin's thickened appearance). This appearance of 
thickening gave credence to the belief that water 
contamination was causing an ice layer to form on the top tube 
between fins. The refrigerant was then changed to R-113 and 
the run was repeated, with the inlet cooling temperature at 
approximately -20 C. Again, a marked thickening of the finned 
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surface appeared on the top of the first tube. The next 
experimental run was conducted on the same finned tube with R- 
113, but at a coolant inlet temperature slightly above 0 C. 
In this situation, no visible thickening of the finned surface 
was observed, but disparities in temperature increases across 
tube lengths persisted. The run was repeated several days 
later to observe whether the phenomenon consistently repeated 
itself. Finally, the smooth copper tube was substituted for 
the finned tube and the experimental run repeated at the 
higher coolant inlet temperature. Noticeable surface 
contamination, having an orange peel appearance, recurred. In 
addition, during this run, there were indications that a 
pocket of non-condensable gases formed in the top of the 
condensing chamber. As every effort had been made to evacuate 
the apparatus prior to commencing the experimental run, and no 
evidence exists to support outside leakage of air into the 
system, the possibility exists that contamination within the 
system produced the non-condensable gases during boiling. 

Two possible sources of the contamination are 
conjectured. The first is that high temperatures on the 
heater tubes and poor circulation in the boiling chamber are 
combining to break down the refrigerant molecules producing 
some type of hydro-carbon. The second is that a chemical 
reaction is occurring between the refrigerants and gasket 
material in the apparatus producing a hydro-carbon. Time and 
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funding limitations have prevented further investigation of 
the problem. 

5 . Summary 

The original test plan called for the testing of 
various enhanced surfaces in a simulated bundle during 
condensation of R-114. The contamination problem described 
above, coupled with the time delays inherent in achieving 
solutions to the other aforementioned encountered problems, 
severely limited the results of this thesis. 

B. DATA REDUCTION 

1. Description of Program Capabilities 

The computer program, DRPIF, that collects and 
processes raw data, in conjunction with a Hewlett-Packard 
computer/data acquisition system described in Chapter II, is 
listed in its entirety in Appendix C. The program is designed 
to calculate and plot heat transfer parameters for a variety 
of different tube bundles, utilizing R-114, R-113, or steam as 
the working fluid. The program has the added capability of 
allowing testing for single tube performance. The program 
consists of five main sections, as follows: 

1. Driver Program, 

2. Main Program, 

3. Property Subroutines, 

4. Modified Wilson Plot Subroutine. 

5. Plotting Subroutines. 
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Of these five sections, only the modified Wilson Plot 
subroutine will be described in lengthy detail. 

The driver program permits the user to take data, 
reprocess data, or plot re-processed data through various 
subroutines. The driver program is listed in lines 1000 
through 1125, in Appendix C. 

The main subroutine (lines 1130-2315) can be divided 
into five parts. The first part allows the user to select the 
physical parameters used in data reduction. The selection of 
parameters consists of the working fluid to be used (R-114, R- 
113, or steam), the vapor saturation temperature to be used 
(derived from averaged thermocouple readings in either the 
vapor section or liquid region of the apparatus) , the 
instrumented test tube type, and whether data will be taken 
on a bundle or individual tubes. The second part allows the 
user to reprocess data with the same selection of physical 
parameters described above, but calls the modified Wilson Plot 
subroutine to calculate the inside and outside coefficients 
used in the correlations. Basic data reduction takes place in 
the third part (lines 1945-2125) . Subroutines are called to 
calculate the properties of the ethylene glycol-water 
solution. These properties, with the physical tube 
dimensions, are used to calculate velocities and Reynolds 
numbers. Low coolant velocities and high viscosities prompted 
the use of twisted tape inserts (thickness 0.559 mm, with a 
pitch for a 180 degree twist of three times the tube's inner 
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diameter) to increase the inside heat flux. The inside 
Nusselt number with a twisted tape insert, provided by Hong 
and Bergles [Ref. 21], is given by the correlation: 



iiS. = ^ = 5.172 



Re 



(4.1) 



where the Reynolds number for coolant flow was given by: 



Re = 4*r/-n'*n* (D. - 46) 
s 1 



(4.2) 



where, 6 is tape thickness. 

The outside heat transfer coefficient is then given by 
the well-known summation of resistances to heat transfer: 



^ + R. + R + — 
U A — , f m — 



(4.3) 
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which algebraically reduces to: 



a = 



1 ^ 
FT 

U mo t o — A- 
o 1 



(4.4) 



The heat transferred to the coolant (Q) is given by the 
relationship: 



Q = 



r-C • (AT) 
p 



(4.5) 
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The heat flux (Q") is subsequently calculated by dividing by 
the outside surface area: 



Q" 



0 



(4.6) 



Where = t^Do ‘L. And, the overall heat transfer coefficient 
(Ug) is given by; 



U 



Q" 

IMTD 



(4.7) 



where the Log Mean Temperature Difference (LMTD) is defined 
to be : 



LMTD = 



AT 



T .-Tc. 
sat 



log [ 



T . -Tc^ 
sat o 



(4.8) 



It should be noted that the wall resistance due to fouling was 
assumed to be negligible for the purposes of calculation. 

The fourth part of the main subroutine creates a raw 
data file (lines 1680-1835) , allowing subsequent reprocessing 
by the modified Wilson subroutine. The fifth part of the main 
subroutine provides a printed output both while taking initial 
data and subsequently after reprocessing data. 

The third major section of the computer program 
calculates fluid properties through called functions for both 
the working fluid and the coolant. The calculated coolant 
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properties as a function of temperature (lines 2330-2615) , are 
kinematic viscosity, specific heat, density, Prandtl number, 
and conductivity. 

The fifth major section of the program provides 
plotting routines for the output files generated in the 
program. The relationships graphically displayed are heat- 
transfer coefficient ratios (either based on the first tube in 
the bundle or as a ratio of the Nusselt value) plotted against 
tube position, heat transfer coefficient plotted against 
either the heat flux or temperature rise of the coolant, and 
the X-Y plot generated from the modified Wilson Plot results. 

2 . Modified Wilson Plot 

The modified Wilson Plot, as outlined by Marto [Ref. 
22], accomplishes an indirect measurement of the outside heat- 
transfer coefficient. The implicit assumption is that the 
overall heat transfer coefficient (Uq) is reliably known from 
data and, therefore a summation of heat transfer resistances 
is assumed. This summation relationship is given by: 



ir ' + Vo + 

o a_ 1 



(4.9) 



where resistance due to wall fouling is assumed equal to zero. 
The summation equation is transformed to a linear 
relationship, as follows: 



[ri -R A ] 
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(4.10) 
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where ; 



a. = C. (^) ; a = C F 

1 1 D. o o 

1 



(4.11) 



which results in the simple linear form, Y = mX + b. It 
should be noted that Theta (defined in line 3290) is derived 
from the Hong and Bergles [Ref. 20] relationship for the 
inside Nusselt number, given by: 



m. = = C. [1 +5.4838-10 ^(Pr®'"^) • (_S)l-25jl/2 

1 A 1 y 



(4.12) 



hence. 



0 = [1 + 5.4838-10 ^(Pr°*'^)-(^)^'^^J^^^ 



(4.13) 



Further, F (defined in line 3330), is derived from Nusselt's 
relationship for a horizontal cylinder subjected to a constant 
heat flux [Ref. 4], given by: 



= .655[A^-p^-g-AhyirD^-Q"]^/^ 



(4.14) 



Hence , 



F = [A^-p^-g-Ahyr)-D^-Q"J^/^ 



(4.15) 



38 



X and Y values are calculated from raw data and the data are 



fit with a least squares approximation. Initial assumed 
values are taken from the aforementioned correlations, and the 
solutions iterated to find the inside coefficient (Ci) and the 
outside coefficient (Cq) that fit the data, where: 



D -F 

Y = )]F ; X = 

U in o A 

o 



( 4 . 16 ) 



the slope (m) is given by: 




and the intercept (b) is given by: 




The accuracy of this method relies heavily upon the number of 
data points taken and the range of velocities utilized. 
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V. RESULTS AND DISCUSSION 



A. EXPERIMENTAL SEQUENCE 

The data runs are summarized in Table 5.1 according to 
tube type, refrigerant used, and approximate operating 
parameters. Specifics for each run are provided in the 
following section under the run title. 



TABLE 5.1 

SEQUENTIAL LISTING OF DATA RUNS 



Run 

Title 


Tube 

Tvoe 


Refrigerant 

Used 


Coolant 
T inlet 




Vapor 

T saturation 


SMT02 


smooth 


R-114 


-20 C 




17.8 


C 


CNFTOl 


finned 


R-114 


-21 C 




18.3 


C 


CNFT02 


finned 


R-113 


-21 C 




64.4 


C 


CNFT03 


finned 


R-113 


-3 to +6 


C 


52.2 


C 


CNFT04 


finned 


R-113 


-3 to +6 


C 


46.3 


C 


SMT03 


smooth 


R-113 


-1 to +6 


C 


56.7 


C 



B. EXPERIMENTAL RESULTS 

1 . Smooth Tubes with R-114 and Low Coolant Temperature 
SMT02 was an experimental run made with four smooth 
copper tubes while condensing R-114. Vapor saturation 
temperature was 17.8 C. Ambient temperature was 2 0 C. The 
ethylene glycol coolant inlet temperature was -20 C. The 
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experiment was conducted following a cleaning of the boiling 
chamber, two flushings of the apparatus, once with acetone and 
once with R-113, evacuation to 27.5 in Hg. , and filling the 
system from an unopened cylinder of R-114. After filling the 
system, the apparatus pressure stabilized at 17 psig. The 
system pressure was lowered to 10 psig by opening the coolant 
flow through the auxiliary condenser. When the desired gage 
pressure was reached, the lowest flow rate was set through the 
instrumented tube bundle, and all three heater units were set 
for power levels corresponding to approximately 10 kw total. 
Gage pressure was maintained nearly constant by controlling 
coolant flow through the auxiliary condenser. The system 
required approximately 35 minutes from lightoff to reach a 
steady state condition, indicated by nearly constant gage 
pressure. The R-114 appeared clear at the commencement of the 
data run. By the time the apparatus reached a steady state 
condition, the contamination matrix was fully formed and 
visible on the top two tubes in the bundle. Data were taken 
at nine different coolant velocities, after the voltage 
indicator on the data acquisition system for the bottom tube 
thermopile appeared to reach a fixed value following each 
velocity change. Data for the tube bundle run (Figure 5.1) 
clearly demonstrates an observed uncertainty of about 7% in 
the overall heat transfer coefficient (Uq) at velocities close 
to 0.9 m/s. This observed uncertainty, in contrast to the 
maximum calculated uncertainty of 4.5% (see Appendix A), 
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Fiaure 5.1 SttT02 — Tube Bundle Performance 



indicates the effect of the contamination on the measurement 
of the overall heat transfer. The most significant result of 
the contamination is the addition of an un-quantif iable 
resistance to heat-transfer that has produced heat-transfer 
coefficients approximately half of those predicted by the 
Nusselt correlation. Further, it is apparent that during 
bundle operation, the contamination degrades the performance 
of the tubes in a graduated manner with the first tube being 
the most adversely affected and the third tube appearing the 
least affected. 

Upon completion of the bundle run, the flow in all of 
the instrumented tubes was shut off. Gage pressure was 
maintained at 10 psig by controlling coolant flow through the 
auxiliary condenser. The coolant inlet temperature remained 
approximately -20 C. The apparatus remained in this condition 
until the tube surfaces appeared to dry off, which took 
normally about ten minutes. Data were then taken at the same 
nine flow rates, as for the bundle, but with flow through only 
one instrumented tube at a time, starting with the top tube. 
Between velocity changes, approximately 10 minutes was allowed 
for the temperature changes in the tube to take effect. Data 
runs were made only after inundation from the previously 
tested tube was no longer present. The R-114 in the boiling 
chamber remained visually clear for all single tube runs. 
Again, as in the bundle data, data from the single tube runs 
(Figure 5.2) fall within an acceptable uncertainty band. 
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Figure 5.2 SMT02 — Individual Tube Performance 



From the conditions under which these runs were made, it is 
noteworthy and logical that the effect of the contamination 
increases in magnitude as a run progresses. The first tube 
gives the lowest results possibly because of residual 
contamination from the bundle run. The first and second tubes 
produced approximately the same results as previously 
discussed in the bundle run, however the third and fourth 
tubes gave lower results when operating as single tubes than 
when operating in a bundle, possibly because condensate 
inundation that provides a rinsing effect in the bundle 
operation, is no longer present during single tube operation. 

2 . Finned Tubes with R-114 and Low Coolant Temperature 
CNFTOl was an experimental run made with copper-nickel 
low integral-fin tubes condensing R-114. Vapor saturation 
temperature was 18.3 C. The ethylene glycol coolant inlet 
temperature was at -21 C. The experiment commenced following 
the completion of the SMT02 single tube runs, transfer of the 
R-114 to the auxiliary storage unit, installation of the 
copper-nickel finned tubes, evacuation of the apparatus to 
27.5 in Hg, and a refill of R-114 from the storage unit. This 
process took approximately six hours to complete. The same 
data taking procedure as outlined for the smooth tube runs was 
followed. The appearance of the R-114 was clear at the 
commencement of the run. By the time the system reached 

steady state, a marked thickening of the fin areas on the top 
tube when compared to the fins of the second tube was apparent 
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(see Figure 4.1). The R-114 appeared clear at the end of the 
bundle run. Data from the bundle run (Figure 5.3) show an 
increase in the overall heat transfer coefficient of 
approximately 100% when compared to the bundle results for the 
smooth tubes SMT02 (Figure 5.1). The coolant velocity 
differences between the smooth tubes and the finned tubes is 
due to the smaller inside diameter of the finned tubes at the 
prescribed mass flow rate. There is no disparity between 
calculated and observed uncertainties given the calculated 
uncertainty band, but performance remains approximately half 
of that expected when compared to the values obtained with the 
Nusselt correlation and enhancement ratios reported by other 
investigations. Once again, the first and second tubes 
demonstrate the greatest degradation in performance due to the 
contamination. The third and fourth tubes either receive less 
contamination or derive enhancement from condensate inundation 
flow. From observations, the unknown contaminant had a higher 
surface tension than R-114, but to what extent the contaminant 
kept fin root areas flooded and subsequently negated the 
enhanced surface effect, remains obscured in the uncertainties 
of the measurements. 

The single tube runs followed the procedures outlined 
for the single tube runs of the smooth copper tubes. However 
at the end of each single tube run, a marked thickening of the 
fin surface at the top of the tube was noticeable. Vapor 
saturation temperature and coolant inlet temperature were the 
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Fiaure 5.3 CNFTOl — Tube Bundle Perf orir>ance 



same as in the bundle run. Data from the single tube runs 
(Figure 5.4) demonstrate an enhancement of approximately 76% 
in the overall heat transfer coefficient when compared to its 
counterpart in SMT02 . The spread in the data is approximately 
7% and, given the calculated uncertainty band, appears 
reasonable. Performance remains lower than expected from 
existing studies, and clearly the magnitude of the 
contamination effect increases as each run proceeded. 

3 . Finned Tubes with R-113 and Low Coolant Temperature 
Upon completion of the single tube runs described in 
CNFTOl, R-114 was transferred to the storage unit. The 
apparatus was drained of residual refrigerant, the system was 
evacuated to 27.5 in Hg, and filled with freshly distilled R- 
113 by using apparatus vacuum to promote flow from the R-113 
container. After filling, system vacuum was at 11 inches Hg. 
The system was brought to a gage pressure of 5 psig 
corresponding to a vapor saturation temperature of 64.4 C. 
Coolant inlet temperature remained the same as described in 
CNFTOl, -21 C. Vapor saturation temperature was monitored 
through the system pressure, and controlled by regulating 
coolant flow through the auxiliary condenser. The same nine 
data points described previously were taken in accordance with 
the stated procedure. The appearance of the R-113 at the 
beginning of the run was clear and clean. Upon the conclusion 
of the run, however, the R-113 had a slight yellowish tinge. 
During the run, the same marked thickening of the fin surface 
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Fiaure 5.4 CNFTOl — Individual Tube Ferformance 



described in CNFTOl, and shown schematically in Figure 4.1 
occurred on the top tube. The data from the bundle run 
(Figure 5.5) demonstrates only a slight improvement in the 
overall heat transfer coefficient when compared to the bundle 
data for SMT02 (Figure 5.1), and in fact represents an 
approximate 80% decrease in the overall heat transfer 
coefficient when compared to the performance of the bundle in 
CNFTOl (Figure 5.3) . In view of the fact that these data were 
collected after transferring the R-114 to the storage unit and 
filling with freshly distilled R-113, it can only be assumed 
that the contamination was freshly created with the R-113 and, 
in fact is associated with the apparatus. It is also 
possible, based on these results, to infer the presence of 
non-condensable gases generated by the contamination, although 
no indications of non-condensable gases were detected during 
the run. In view of these facts, the reliability of the 
calculated values for the overall heat transfer coefficient is 
extremely doubtful even though data spread falls within the 
acceptable uncertainty band. 

Upon completion of the bundle run, data were taken for 
each of the tubes individually at the same system conditions 
described for the bundle run, and in accordance with the 
procedure described for single tube runs. The appearance of 
the R-113 at the completion of the single tube runs was 
slightly-darkened, with a yellow tinge, but otherwise appeared 
free from any floating contaminants. On each of the single 
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Fiaure 5.5 CNFT02 — Tube Bundle Performance 



tube runs, a marked thickening of the fin surface developed in 
the course of the run. There exists reasonably good agreement 
between the single tube data (Figure 5.6) and the single tube 
data from CNFTOl in the calculation of the overall heat 
transfer coefficient, and no disparity exists between the data 
spread and the calculated uncertainty band. The apparatus was 
evacuated for five minutes prior to beginning the single tube 
runs and if non-condensable gases were present during the 
bundle run, this explains the improved performance of the 
tubes during individual runs. 

4 . Finned Tubes with R-113 and High Coolant Temperature 
This data run was made with the same copper-nickel 
tubes described in CNFTOl and CNFT02 . The system was allowed 
to sit overnight at 11 in Hg, while the sump temperature was 
brought to an operating range of -3 C to +6 C. No change in 
the system vacuum over this 12 hour period occurred. The 
decision to run the system at a higher coolant temperature was 
motivated by a belief that the contamination of the system was 
water, and that at the colder coolant temperature, an ice 
sheath was forming on the tubes inhibiting heat transfer. The 
appearance of the R-113 before commencement of the bundle run 
was a medium yellow tinge. The system was brought to a steady 
state pressure corresponding to a vapor saturation temperature 
of 52.2 C. Difficulty was encountered in maintaining positive 
system pressure with the auxiliary condenser operating. The 
auxiliary condenser was therefore turned off and steady state 
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was maintained by slight power adjustments to the heating 
elements in the boiling chamber. Nine sets of data were taken 
in accordance with the procedure for bundle runs previously 
described. No observable thickening of fin surfaces during 
the run occurred. The data from the bundle run (Figure 5.7) 
shows a slight increase in overall heat transfer coefficient 
over the bundle results obtained in CNFT02 (Figure 5.5) and 
this is possibly explained by a decrease in the viscosity of 
the contaminant when exposed to the warmer tube surface. 

Single tube runs were made on each of the tubes, in 
accordance with procedures outlined previously, at the same 
coolant inlet temperature as in the bundle run, but at a 
system pressure corresponding to a vapor saturation 
temperature of 64.4 C. No change in the appearance of the R- 
113 was detected. Again no thickening of the fin surfaces was 
detected. The single tube results (Figure 5.8) demonstrate a 
corresponding range of values in the overall heat transfer 
coefficient when compared to the results for single tubes in 
CNFT02 (Figure 5.6), but is significantly higher than the 
bundle results. It is possible that non-condensable gases 
were present in bundle operation, but were eliminated by 
evacuation prior to single tube data being taken. No 
indications of non-condensable gases were detected during 
bundle or single tube operation. 
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Fiaure 5.7 — CNFT03 — Tube Bundle Performance 
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VEL. (M/s) 

Fiaure 5.8 — CNFT03 — Individual Tube Performance 



5 . Finned Tubes with R-113 and High Coolant Temperature 

This data run was made in order to see if the results 
produced in CNFT03 were repeatable. The system operating 
pressure corresponded to a vapor saturation temperature of 
46.3 C, and the coolant inlet temperature remained unchanged 
from the coolant temperature in CNFT03 . Nine data sets were 
taken in accordance with the bundle run procedures outlined 
previously. No thickening of fin surfaces during the run was 
observed. The color of the R-113, at the completion of the 
run, was amber. The bundle data (Figure 5.9) does demonstrate 
that the results of CNFT03 were repeatable at least within the 
calculated uncertainty. 

Single tube runs were made in accordance with the 
procedures discussed previously. Nine data sets were taken on 
each tube, at the same conditions specified in the bundle run. 
No thickening of fin surfaces was observed during the runs, 
and the R-113 color remained amber. There is no major 
disparity between the single tube results (Figure 5.10), and 
the single tube results for CNFT03, given the calculated 
uncertainty band. The single tube results were significantly 
higher than bundle performance, and the build up of non- 
condensable gases during bundle operation and evacuation of 
the apparatus prior to commencing single tube runs is 
conjectured as a possible explanation. At no time during 
bundle or single tube operation were indications of non- 
condensable gases detected. 
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Figure 5.9 — CNFT04 — Tube Bundle Performance 
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6 . Smooth Tubes with R-113 and High Coolant Temperature 
SMT03 was a data run made on the same aforementioned 
copper tubes, but with R-113. Vapor saturation temperature 
was 56.7 C. Coolant inlet temperature ranged between -1 C and 
+6 C. The experiment was conducted following a test run with 
R-113 and finned tubes and subsequent replacement of the 
finned tubes, and evacuation to 11 in Hg. While the tube 
change was being made, the R-113 remained in the boiling 
chamber exposed to the atmosphere for a two hour period. The 
R-113 appeared light amber in color at the commencement of the 
bundle run. The desired gage pressure was maintained by 
controlling the boiling chamber heaters through their 
corresponding variac controllers. The auxiliary condenser was 
not utilized because the additional coolant flow provided by 
the auxiliary condenser made it difficult to maintain the 
desired system pressure. During the data run, the top two 
instrumented tubes assumed an orange peel textured or dimpled 
appearance. The R-113 had a darker amber appearance at the 
completion of the bundle run. Again, data were taken at the 
same nine coolant flow settings following the procedure 
outlined previously. Data taken during bundle operation 
(Figure 5.11) , while falling within the acceptable uncertainty 
bands, clearly demonstrates the presence of non-condensable 
gases as well the presence of the contamination acting to 
inhibit the overall heat-transfer performance. The top tube's 
very poor performance can be attributable to its location in 
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Fiaure 5.11 — SI1T03 — Tube Bundle Performance 



a pocket of non-condensable gases. And the bundle performance 
as a whole, when taken in comparison with the performance 
demonstrated in the bundle run of SMT02, indicates an 
increased concentration of the contamination and the presence 
of non-condensable gases. 

Upon completion of the bundle run, coolant flow 
through the instrumented tubes was shut off and gage pressure 
maintained by controlling flow through the auxiliary 
condenser. While no coolant flowed through the instrumented 
tubes, the textured appearance noted during the bundle 
operation gradually transformed to small droplets on the top 
two tubes which then coalesced into larger ellipsoid shaped 
droplets that were held on the tube by surface tension forces. 
Data was taken on the first two tubes in the bundle following 
the procedure outlined previously for single tube runs. At 
the completion of the second run, it was noticed that the top 
view port flanges were cold to the touch while the bottom view 
port flanges were warm to the touch. The assumption was 
therefore made that non-condensable gases had collected in the 
top of the condensing chamber. Before evacuation of the 
apparatus could be effected, the pyrex glass in the top center 
view port cracked. As a result, the apparatus was secured and 
no further data taken. It should be noted that, prior to the 
view port breaking, there was no evidence to suggest that the 
apparatus was other than air tight, so that the source of non- 
condensable gases was most probably due to their generation 
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within the boiling chamber. Data taken in the single tube 
runs (Figure 5.12) reflects not only the effects of the 
contamination, but also the effects of the presence of non- 
condensable gases as well. This data showed the greatest 
variance between calculated and observed uncertainties in the 
overall heat transfer coefficient. This data is limited to 
the performance of only the top two tubes, due to the afore- 
mentioned breakage of the view port and the presence of non- 
condensable gases. 

C. SUMMARY OF RESULTS 

The comparison between finned and smooth tubes shows an 
enhancement ratio in the overall heat transfer coefficient for 
the finned tubes compared to the smooth tubes of approximately 
2.0. While this enhancement in the overall heat-transfer 
coefficient might translate to a vapor-side heat-transfer 
coefficient enhancement ratio of approximately 4.0 to 4.5, 
based upon the calculated values for the vapor-side heat- 
transfer coefficient for the tested smooth tubes when compared 
against the Nusselt's correlation for smooth horizontal tubes 
and the reported enhancement ratios from other finned-tube 
investigations [Refs. 5-13], it is estimated that the effect 
of the contamination during these tests has been a degradation 
of up to 50% in the vapor-side heat-transfer coefficient. The 
contamination has resulted in the promotion of another 
resistance to heat transfer that has no predictable 
characteristics and can not be quantified. This contamination 
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Fiaure 5.12 — SMT03 — Individual Tube Performance 



must therefore be removed before any quantifiable heat- 
transfer results can be obtained. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Based on the data gathered during this investigation for 
condensation of refrigerants R-114 and R-113 in the multi-tube 
apparatus, the following conclusions are reached: 

1. The construction of a multi-tube condensation test 
apparatus begun by Murphy [Ref. 3] and Zebrowski [Ref. 
4] has been completed, instrumented and integrated with 
the various support systems. 

2. Data reduction software has been developed that produces 
high confidence results from experimental readings. 

3. Data were collected from condensation experiments with 
R-114 and R-113 on smooth and finned tubes demonstrating 
an enhancement ratio of approximately 2.0 in the overall 
heat transfer coefficient for the finned tubes over the 
smooth tubes. 

4 . Serious apparatus internal contamination has degradated 
the value of the overall heat-transfer beyond acceptable 
limits, making further data reduction meaningless. 



B. RECOMMENDATIONS 

Based upon the data gathered during the operation of the 
multi-tube condensation test apparatus, and the experience 
gained in construction of this apparatus, the following 
recommendations are made; 

1. Samples should be taken from the R-114 and R-113 used 
in the experimental runs and a chemical analysis 
performed on these samples to assist in determining the 
source of the contamination in the system. 

2. Modifications to the apparatus, as dictated by the 
results of the chemical analysis of the refrigerants, 
should be made. 
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3 . 



New flow meters allowing a greater range in coolant 
velocities should be installed. The greater range in 
coolant velocities will enhance the accuracy of the 
modified Wilson plot method in calculating the outside 
heat-transfer coefficient. 

4 . The temperature control system for the eight ton 
refrigeration unit should be re-evaluated to determine 
if tighter control of sump temperature can be 
maintained . 

5. After appropriate modifications to the apparatus have 
been made, the smooth tube and finned tube experiments 
should be repeated and contrasted with the results 
presented here, in order that baseline tube perfomnance 
might be established. 

6. During operation of the system, flow should be throttled 
to the coolant supply pumps with the valves on the 
coolant inlet side. This will increase the pump life 
and reduce noise in the laboratory space. 
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APPENDIX A 



UNCERTAINTY ANALYSIS 



A certain amount of uncertainty exists in any engineering 
measurement. These uncertainties arise both from known 
sources, such as calibration and measurement errors from 
sensing devices, and also unknown sources, such as operator 
experience and unexpected experimental reactions. The 
procedure used in this uncertainty analysis is based on the 
Kline-McClintock [Ref. 21] method. This method assumes that 
a result R is a function of the variables that contribute to 
that result. Therefore the uncertainty in a result is a 
function of the uncertainties in each of the variables. This 
method is expressed mathematically by: 






(A.l) 



The measurement of the overall heat transfer coefficient 
(Uq) performed in this thesis is given by: 



^o A .I2OTD 
o 

where the heat transferred (Q) is given by: 



(A. 2) 



Q = r-C • (AT) 
p 



(A. 3) 
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and the Log Mean Temperature Difference (LMTD) is defined as: 



IKTD = 



(AT) 



T . -Tc . 

, r sat 1 , 
log[= — 



T . -Tc 
sat o 



(A. 4) 



Therefore the uncertainty in the measurement of the overall 
heat transfer coefficient (Uq) is given by: 



6U 6C o xA-n. o <5D _ 

° = [ (^) 2 + 2 ^ 2 2 



U 



AT' 



-6LMTD> 2 1/2 
^ IMTD^ ^ 



(A. 5) 



6A 

where — 5. 
A 
o 

given by: 



reduces 




and the uncertainty in the LMTD is 



6IMTD 

LMTD 



lA^ + + 






(A. 6) 



where A is given by the expression: 



A = ST^ [ 



AT 



satMT^^-Tc.) 



log [ 



T . -Tc . 
sat i- | 

T . -Tc ■' 
sat o 



] 



where B is given by the expression; 



B = 



6Tc^ 






T ,-Tc. 
sat 



log I 



T . -Tc 
sat o 



where C is given by the expression; 
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1 



c = 



6Tc 

o 

(T . -Tc ) 
sat o 



log [ 



T . -Tc . 
sat 1, 

T ^-Tc ^ 
sat o 



It should be noted that the uncertainties for mass flow 
rate (p) and the coolant temperature rise across the 
condensing tubes (At) were calculated during calibration of 
the flow meters and thermopiles respectfully. Reference 
material providing data on thermophysical properties, 
specifically the specific heat (Cp) , listed no uncertainties 
associated with the curves for aqueous-ethylene glycol 
solutions. The composition of the aqueous-ethylene glycol 
solution by weight percentages was calculated by measurement 
of the specific gravity and the change in this composition 
over the experimental period was determined to be negligable 
and therfore this uncertainty has been ignored. 

The uncertainty in the overall heat transfer coefficient 
represents the maximum uncertainty present from the 
calibration of the four flow meters and the maximum 
uncertainty in the calibration of the four thermopiles, 
regardless of association, and therefore represents a 
conservative estimate of the uncertainty present in each 
calculation of the overall heat transfer coefficient (U^) . 
This maximum uncertainty was calculated to be 4.5 percent. 

The disparity between the calculated uncertainties and the 
uncertainties observed by the spread in data particularly in 
the single tube runs, can be accounted for by the as-yet 
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undetermined chemical or phase reactions generated by the 
contamination in the system. 

The values of the uncertainties in the measured variables 
that were used to calculate the uncertainty by the Kline- 
McClintock method are listed in Table A.l. 



VARIABLE 

8T_ 

r 

5C 
P 

5 

6 AT 
AT 



TABLE A.l 

UNCERTAINTY VARIABLES 

VALUE REMARKS 

0.01 From flow meter 

calibration data 



0.00 Not available 



0.01 From thermopile 

calibration data 



5D 



o 



D 



o 



0.002 From micrometer 

name plate data 



6LMTD 

LMTD 

SMT02 

CNFTOl 

CNFT02 

CNFT03 

CNFT04 

SMT03 



Listed below by 
run title 

0.005 

0.002 

0.001 

0.003 

0.002 

0.017 
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APPENDIX B 



LIGHTOFF AND SECURING PROCEDURES 
A. SYSTEM LIGHTOFF 

1. Push the starter button in the control box for the re- 
circulation pump. This control box is located on the bulkhead 
above the re-circulation pump in the outside area adjacent to 
the refrigeration unit. 

2 . Turn the switch on the refrigeration unit control 
panel, located in front of the refrigeration unit to the 
"auto" position after passing through "on" position. 

3. Set the desired temperature on the roughly graduated 
Fahrenheit scale on the control panel thermostat. It requires 
approximately four hours to chill the sump 40 degrees C. The 
thermometer located on the side of the ethylene glycol sump 
must be monitored to ensure the desired sump temperature is 
attained and maintained. Slight adjustments in the 
refrigeration unit thermostat can be expected due to the 
coarseness of its scale. 

4 . Energize the heater variacs by switching on the 
breakers in the breaker panel located in the laboratory space 
on the bulkhead next to the counter. 

5. Set the heater variacs to the desired position, after 
ensuring that the switch panel for the heater tubes located on 
the apparatus has all switches in the "on" position. Monitor 
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apparatus pressure through the pressure gage at the top of the 
apparatus, ensuring system pressure does not exceed 30 psig. 

6. Turn on the pump motors by pushing down on the arm of 

the appropriate breaker box for the pumps located on the 
bulkhead next to the ethylene glycol sump. The pumps are 
marked ••auxiliary condenser^^ and ••instrumented tube 

condenser, •• respectively. Flow in the auxiliary condensate 
system can be controlled with the individual gate valves 
located at the coil penetrations on the apparatus. The 
auxiliary condenser will produce the fastest adjustments to 
system pressure if pressure is rising too quickly. Flow 
through the instrumented tubes can be controlled by the ball 
valves located at the bottom of the respective flow meter. 

7 . Throttle down coolant flow to the supply pumps with 
the valve located before the pump suction. This will increase 
the pump life and reduce noise in the laboratory space. 

B. SECURING PROCEDURES 

1. Turn all variacs to the zero position and switch off 
all breakers in the power panel on the bulkhead. 

2. Turn the breakers for the pumps to the off position at 
the switch boxes near the ethylene glycol sump. 

3. If apparatus will not be operated for an extended 
period, turn the switch on the refrigeration control panel to 
the ••off^^ position after passing through ••on.^^ 

4. Allow the re-circulation pump to operate for at least 
ten minutes after switching off the refrigeration unit to 
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dissipate any back pressure in the system; then secure the 
pump. 
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APPENDIX C 



1000 I 
1005 » 
1010 I 
1015 » 
1020 ' 
1025 I 
1030 
1035 
1040 
1045 
1050 
1055 
1060 
1065 
1070 
1075 
1030 
1085 
1090 
1095 
1100 
1105 
1110 
1115 
1120 
1125 
1130 
1135 
1140 
1145 
1150 
1155 
1160 
1165 
1170 
1175 
1180 
1185 
1190 
1195 
1200 
1205 
1210 
1215 
1220 
1225 
1230 
1235 
1240 



DATA REDUCTION PROGRAM 



FILE: 

PURPOSE: 



DRPIF 

This program collects and processes condensation data fo'' 
the R-114 tube-bundle apparatus. 

NOUEMBER 2, 1988 



SELECT OPTION 

0 TAKING DATA OR REPROCESING PREVIOUS DATA“ 

PLOTTING H VS DELTA-T 

PLOTTING HRAT VS N""" 

PLOTTING UILSON 

PURGE FILES- 
XYREAD 

NUSSELT ESTIMATE 



CREATED 
UPDATED: November 8, 1988 
CHANGE H0D=3 AND USING CALCULATED Cl FROM UILSON 
BEEP 

PRINTER IS 1 
PRINT USING "4X, 

PRINT USING "BX, 

PRINT USING “6X, “”1 
PRINT USING ''BX/'"2 
PRINT USING “EX, *‘•'3 
PRINT USING ”BX , “ ‘’4 
PRINT USING "GX, *•"5 
PRINT USING "6X/'”6 
PRINTER IS 701 
INPUT Icall 

IF Icall-0 THEN CALL Main 
IF Icall-1 THEN CALL Plot2 
IF Icall-2 THEN CALL Plotl 
IF Icall=3 THEN CALL Plot3 
IF Icall-4 THEN CALL Purge 
IF Icall=5 THEN CALL Xyread 
IF Icall*=6 THEN CALL Nusselt 
END 

SUB Main 
COM /Cc/ C(7) 

COM /Fid/ Ift 

COM /Nus/ 1 1 n , Tsai , Odp 1 , Hnus 

COM /Uill/ Doa(4 ),Dia(4 ),Kma(4 ), lact 

COM /Uil2/ Del t a , I sat , Nse t s , Hod , Cia( 3 ) , Alpaa( 3 ) 



DIM Fma(8,4 ),Fm(4),Emf(8),Tp(3),T(8 ),Ho(3),0dp(3),Uo(3) 

DATA 33.4,20,34.4,34.5,23.5 

DATA 41.3,25,41.3,41.6,28.8 

DATA 49.1,30,48.4,48.6,34.1 

DATA 57.0,35,55.3,55.6,38.4 

DATA 64.8,40,62.3,62.7,44.6 

DATA 72.7,45,69.2,69.7,48.8 

DATA 80.6,50,76.3,76.7,55.2 

DATA 88.4,55,83.3,83.8,60.5 

DATA 96.3,60,90.2,90.8,65.6 

READ Fma( * ) 

DATA 5.172,5.172,5.172,5.172 
READ Cia( ♦ ) 

DATA 0.10086091,25727.94369,-767345.8295,78025595.81 
DATA -9247486589,6. 97B88E11 ,-2.661 92E13, 3. 94078E1 4 
READ C( • ) 

DATA 0.015875,0.014000,0.0,0.0,0.0 
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125C3 

1255 

1250 

1255 

1270 

1275 

12B0 

1255 

1290 

1295 

1300 

1305 

1310 

1515 

1320 

1325 

1330 

1335 

1340 

1345 

1350 

1355 

1360 

1365 

1370 

1375 

1360 

1355 

1390 

1395 

1400 

1405 

1410 

1415 

1420 

1425 

1430 

1435 

1440 

1445 

1450 

1455 

1460 

1465 

1470 

1475 

14B0 

1485 

1490 

1495 

1500 

1505 

1510 

1515 

1520 

1525 

1530 



DATA 386.0,42.975,0.0,0.0,0.0 
DATA 0.0005598,3 

READ Doa( * ) , D 1 a ( ♦ ) , Kna 1 * ) , De 1 1 6 , Hoc ' HoC=H/Di 
L=1.2192 ! Condensing length 

Jset *0 
BEEP 

INPUT "ENTER MONTH, DATE AND TIME ( MM DD • Hh ■ MM • S5 )", Dt gS 
OUTPUT 709i "TD" ;Dtg$ 

BEEP 

INPUT "SELECT OPTION (0 = DAQ, l=FlLE)'‘,In 

I hard= 1 

BEEP 

INPUT "WANT A HARDCOPY PRINTOUT ( 1 =DEF = YES , 0-N ) " , I hard 
BEEP 

INPUT "SELECT ( 0-R- 1 1 4 , 1 =STEAM , 2 =R- 1 1 3 , 3 = EG ) " , I f t 

I in=l 

Isat-2 

BEEF 

INPUT "SELECT SAT TEMP MODE ( 0 = L I Q , 1 =UAF . 2 = ( L I 0-^UAF )/2=DEF ) " , I sat 
IF Ihard'l THEN PRINTER IS 701 
IF In-0 THEN 
BEEP 

INPUT "GIUE A NAME FOR THE NEW DATA FILE",File$ 

CREATE BOAT File$,20 
ASSIGN @F:le TO FileS 
BEEP 

INPUT "ENTER TUBE CODE",Itube 
BEEP 

INPUT "ENTER EG CONCENTRATION ( WT PERCENT )", Egra t 
ENTER 709iDtg$ 

OUTPUT @File,Dtg$ 

OUTPUT @File; Itube,Egrat ,Ddl ,Dd2,Dd3.Dd4,Dd5 

I act*=10 

BEEF 

INPUT "SELECT ( 0=TOP, 1=SEC0ND 1 0 = BUNDLE = DEFAULT ) ",I ac t 

PRINT 

PRINT USING "10X,“"f^ILE NAME'- " " , 1 2 A " . F 1 1 e$ 

PRINT 

ELSE 

BEEP 

INPUT "ENTER NAME OF EXISTING FILE' ,F:le$ 

BEEP 

INPUT "ENTER NUMBER OF DATA SETS STORED ", Nse t s 

I wi 1*1 

BEEP 

INPUT "WANT TO CALL WILSON? (1 =DEFAULT* YES , 0*NO ) " , I u 1 1 
IF Iwil*l THEN 
BEEP 

INPUT "WHICH TUBE ( 0=TOP , 1 =SECOND , . . , 10 = BUNDLE ) " , I ac t 
CALL Wilson 

IF Ihard*l THEN PRINTER IS 701 
PRINT 

PRINT USING “10X,""FILE NAME: " " , 1 2A" ; F 1 1 e$ 

PRINT 

IF Iact>9 THEN 

PRINT USING "10X, ""INSIDE COEFFICIENTS: " " , 4 ( DO . 3D , 2X ) “ . Ci a( O 
PRINT USING "10X, ""ALPHAS' " " , 4 (1 X , M2 . 3DE , 2 x ) " . A 1 p 
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aa( ♦ ) 

1535 

1540 

0" i lac 

1545 

D " ; lac 

1550 

1555 

1565 

1570 

1575 

1580 

1585 

1586 
1590 
1595 
1600 
1605 
1610 
1615 
1620 
1625 
1630 
1635 
1640 I 
1645 
1650 
1655 
1660 
1665 
1670 
1675 » 
1630 
1685 
1690 
1695 
1700 
1705 
1710 
1715 
1720 
1725 
1730 
1735 
1740 
1745 
1750 
1755 
1760 
1765 
1770 
1775 
1780 
1785 
1790 
1795 
1800 
1805 I 

. 1810 



ELSE 

PRINT USING "10X, "“INSIOE COEFFICIENT FOR TUBE "\0,“” 

t , Cia( lac t ) 

PRINT USING “10X , “"ALPHA FOR TUBE "",0/”‘^ 
t , Alpaa( lact ) 

ENO IF 
ELSE 

INPUT "ENTER Cl UALUES ( 0EFAULT=5 . 172 )" , Cia( * ) 

END IF 

ASSIGN ©File TO FileS 
ENTER ©File;Otg$ 

ENTER ©Fileiltube, Egrat , Ddl , Od2 , Dd3 , Dd4 , Dd5 
INPUT “ENTER TUBE COOE’Mtube 
END IF 
Iout*l 
BEEP 

INPUT "UANT TO CREATE AN OUTPUT FILE7 ( 1 =DEF=YES , 0=N0 ) ",Iout 
IF Iout=l THEN 
BEEP 

INPUT "ENTER A NAME FOR OUTPUT FILE",Fout$ 

CREATE BOAT Fout$,5 
ASSIGN ©Fout TO FoutS 
END IF 



Oo*Doa( I tube ) 

Oi*Oia( I tube ) 

Kn=Kr)a( 1 1 ube ) 

Ax*PI*0i"2/4 » Cross-sectional area 

Ao*PI ♦0o*L 

Rr)«Do*L0G( Oo/Di )/(2*Kn) 



IF ln*0 THEN 
PRINTER IS 1 
BEEP 

PRINT "SET FLOWMETER REAOINGS CORRESPONOING TO:" 

PRINT " " iFna( Jset , 1); "7. OF METER 2 ANO HIT CONTINUE" 

PAUSE 

PRINTER IS 701 

OUTPUT 709; "AR A«^0 ALB UR5" 

Nend=8 » INCREASE TO 9 IF FIUE TUBES IN BUNOLE 
FOR 1=0 TO Nend 

OUTPUT 709; "AS SA" 

Usun=0 

FOR J=1 TO 5 
ENTER 709; E 
Usun=Usun+E 
NEXT J 

Enf ( I >=Usun/5 
NEXT I 

OUTPUT 709; "AR AF20 AL23 UR5" 

FOR 1=0 TO 3 

OUTPUT 709; "AS SA" 

Usun*0 

FOR J=1 TO 50 
ENTER 709;E 
Usum=Usun+E 
WAIT .25 
NEXT J 



,00.3 

,00.3 
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181E 

1620 

1825 

1830 

1835 

1840 

1845 

1850 

1855 

1860 

1885 

1870 

1875 

1880 

1890 

1895 

1896 

1900 

1905 

1906 

1907 

1906 

1910 

1915 

1920 

1925 

1930 

1935 

1940 

1945 

1950 

1955 

1960 

1965 

1970 

1975 

1980 

1985 

1990 

1995 

2000 

2005 

2010 

2015 

2020 

2025 

2030 

2035 

2040 

2045 

2050 

2055 

2060 

2065 

2070 

2075 

2080 

2085 



Tp ( I ) = 'J5un/50 
NEXT I 
ELSE 

ENTER @Fiie.Fn( * ),Enf ( ♦ ),Tp( ♦ ) 

END IF 

DATA ANALYSIS 
Nend=8 

FOR 1*0 TO Nend 

T( I )*FNTvsv(Enf ( I ) ) 

NEXT I 

Tvap-(T(0 ) + T( 1 ) + T(2 ) )/3 
.Tliq*(T(3) + T(4) )/2 
T sat*T 1 iq 
END IF 

IF Isat'l THEN 
T sat'Tvap 
END IF 

IF Isat=2 THEN 

Tsat*(Tvap+TI:q)/2 
END IF 

FOR 1*0 TO 4 

Fn(I )*Fma(J5et,I ) 

NEXT I 
Jset* Jset+1 
PRIN"^ 

PRINT USING "10X,""Data set number * "^DD'^iset 

PRINT 

lbeg*0 

Iend=3 

IF Iact'10 THEN 
I beg* I ac t 
lend* I ac t 
END IF 

FOR I*Ibeg TO lend 

Grad = FNGrad( Em f ( 1+5 ) ) 

Delt*ABS(Tp( I )/(Grad*10 ) ) 

Tavg=T ( 1+5 )+Deit* .5 
Rhoeg*FNRhoeg( Ta vQ , Egra t ) 

Nueg*FNNueg( Tavg , Egrat ) 

Mueg=Nueg*Rhoeg 
Cpeg*FNCpeg( Tavg , Egrat ) 

Keg*FNKeg( Tavg , Egrat ) 

Preg*Cpeg*Hueg/'Keg 
Mdot*FNFmcal( I ,T( 1-^5 ),Fmi I ) ) 

Ueg=Mdot / ( Rhoeg* A>, ) 

Reeg*Ueg*Di/Nueg 

Re5 = 4*Mdot/(PI*f1ueg*(Di-4*Deita ) ) 
Qdot*Mdot*Cpeg*Deit 
Odp( I )=0dot /Ao 
IF 1*0 OR I=Iact THEN 
Odp 1 =0dp ( I ) 

CALL Nusseit 
END IF 

Lmtd*Delt/LOG( (Tsat-T(I+5) )/ ( Tsat-T( I +5 )-Dei t ) ) 
Uo( I )*0dp( I )/Lmtd 
IF Reeg<4000 THEN 
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2090 

2095 

2100 

2105 

2110 

2115 

2120 

2125 

2130 

2135 

2136 
) 

2137 
2140 

2145 

2146 

2147 

2148 

2150 

2151 

2152 

2155 

2160 

2161 

2165 

2166 
2167 
2170 
2175 

HO 

2180 

K ) 

2185 

,3X,r 

2190 

2195 

,3X" 

2200 

2205 

2210 

2215 

2220 

2225 

2230 

2235 

2240 

2245 

2250 

2255 

2260 

2265 

2270 

2275 

2280 

2285 

2290 

2295 

2300 

2305 



Nueg = C 13' I )♦( l+5.484E“3*Preg".7»<Res/Hod)''1.25)''.5 
ELSE 
BEEP 

PRINT USING " 10X INCORRECT TURBULENT CORRELATION""" 
Nueg-.027*Reeg^ ,8*Preg" .3333*Cfeg 
END IF 

Hi*Nueg*Keg/Di 

Ho(I )=!/( 1/Uo^ I )-Do/(Di*Hi )-Rn ) 

IF 1=0 OR I=Iact THEN 

PRINT USING "10X,""Mass flow rate “ " MZ . 3DE " ; Mdot 

PRINT USING " 10X Inside Tube Dia. (n.) = " “ , MZ . 3DE " i 0 i a( I t ube 



PRINT USING "10X, 
PRINT USING "10X, 
PRINT USING ”10X, 
PRINT USING "10X, 
PRINT USING "10X, 
PRINT USING "10X, 
PRINT USING "10X, 
PRINT USING "10X, 
PRINT USING "10X, 



"180 DEG OUER Dia. (HOD) 
"Inlet tenperature 
"Saturation tenp (Deg C) 
"DELTA Tape Thickness 
"DELT temp Dif. 

"Log. Mean tenp Dif. 
"Heat flux 
"Conductivity E.G. 
"Conductivity Tube Metal 



•• ,MZ.3DE" ;Hod 
" " .MZ.3DE" ;T( 1+5 ) 
" " ,MZ .3DE" ;Tsat 
"",MZ.3DE"i Delta 
" " ,MZ.3DE" ;Delt 
"",MZ.3DE";Lntd 
"".MZ.3DE";Qdp(I ) 
"",MZ.3DE";Keg 
"" ,MZ.3DE" ;Kn 



PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

HNUS 

PRINT 



USING 

USING 

USING 

USING 

USING 

USING 



■10X/ 
"10X, ' 
•10X, ’ 
"10X, ' 
•10X, • 
'■ 10X , ' 



USING "10X, 
USING "10X, 



"Prandtl nunber 
"Reynolds numcer 
"Reynolds number H&B S 
"Inside h , t . c . 

"Inside NUSSULT NO. 
"OUERALL H.t .c. ( Uo ) 



"TUBE 

t 



"",MZ.3CE";Preg 
" " , MZ . 3DE " , Reeg 
" " ,MZ.3DE" iRes 
" ‘\MZ.3DE" iHi 
" " , MZ . 3DE " i Nueg 
"" ,MZ.3DE" iUo( I ) 



ELSE 



FM 


UEG 


DELT 


Uo 




( X ) 


( n/ 5 ) 


( K ) 


( U /m 




Lg 

X 


.DD,4X,DD. 


.3D,2X,M2. 


.3DE,3X,M2. 


3DE 


I ) , Hnus 








I,3X,Z 


.DD.4X,DD, 


.3D,2X,MZ, 


.3CE,3X ,m: . 


3DE 



I + l ,Fm( I ) ,Ueg,Deit ,Uo( I ) ,Hc( I ) 
END IF 
NEXT I 



IF lm=0 THEN 



Okaccpt=l 

BEEF 

INPUT "OK TO ACCEPT THIS SET ( 1 =DEFAULT = YES , 0 = NO , Okaccp t 
IF Okaccpt = l THEN OUTPUT @F 1 1 e . Fm (♦ ) , Enf ( * ) , Tp ( ♦ ) 

END IF 

IF (Okaccpt=l OR Im=l) AND Iout=l THEN 
FOR 1=0 TO 3 

OUTPUT @Fout .Hod ),Qdp( I > 

NEXT I 
END IF 

IF lm=0 THEN 
Okrpt= 1 
BEEP 

INPUT "WILL THERE BE ANOTHER RUN ( 1 = YES=DEFAULT , 0 = NO ) " , Okrpt 
IF Okrpt=l THEN 1680 
ELSE 

IF Jset<Nsets THEN 1680 
END IF 

ASSIGN ©File TO ♦ 
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r3i0 

2315 

233C 

2325 

2330 

2335 

2340 

2345 

2350 

2355 

2360 

2365 

2370 

2375 

2380 

2385 

2390 

2395 

2400 

2405 

2410 

2415 

2420 

2425 

2430 

2435 

2440 

2450 

2455 

2460 

2465 

2470 

2475 

2480 

2485 

2490 

2495 

2500 

2505 

2510 

2515 

2520 

2525 

2530 

2535 

2540 

2545 

2550 

2555 

2560 

2565 

2570 

2575 

2580 

2585 

2590 

2595 

2600 

2605 



IF Iout=l THEN ASSIGN @File TO ♦ 

SUBEND 

DEF FNGrad(T) 

Grad=-3.877857E-5-2*4.7142857E-8*T 

RETURN Grad 

FNEND 

DEF FNKcu(T ) 

OFHC COPPER 250 TO 300 K 
Tk=T+273. 15 
K = 434-. 112*Tk 
RETURN K 
FNEND 

DEF FNNueg(Tc,Egr ) 

RANGE OF VALIDITY: -20 TO 20 DEG C 

Tk-Tc+273.15 

Nul =7. 11 96507E-3-Tk*( 7. 4863347E-5-TkM 2.62 94943E-7-Tk *3.08333296-10 ) ) 

Nu2>"4.9237638E-3-Tk*( 4.92l3912E-5-Tk*( 1 . 6437534E-7-Tk*l . 8333331E-10 ) > 

Nu3*8. 6586293E-3-Tk*( 8. 8837902E-5-Tk*( 3. 0495032E-7-Tk *3.49999966- 10 ) ) 

A2 = (Nu3-2*Nu2 + Nul )/200 

Al»(Nu2-Nul-940*A2 )/10 

A0=Nul-42*Al-1764*A2 

Nu*A0 + Egr*( Al+Egr*A2 ) 

RETURN Nu 
FNEND 

DEF FNCpeQ(Tc.Egr) 

RANGE OF VALIDITY: 0 TO 20 DEG C 
Tk-Tc+273. 15 

Cp2*l .4748125E + 3 + TK *6.25 
Cp3*9.5800500E42+Tk*7.3 
A2-(Cp3-2*Cp2 + Cpl )/200 
Al = (Cp2-Cpl-900*A2 )/10 
A0«Cpl-40*Al-l600*A2 
Cp*A0 + Egr ♦ ^ A 1 +Egr * A2 ) 

RETURN Cp 
FNEND 

DEF FNRhoeg( T . Egr ) 

Rol*l .0607093E + 3-TM 3. 7051 283E-1+T *4.08371836-3 ) 

Ro2 = l .0748272E-^3-T*( 4 . 4266 1 95E- 1-^T*4 . 0939706E-3 ) 

Ro3-l .0885934E+3-T*(5.7355653E-l + T*6 . 1281405E-3 ) 

A2 = (Ro3-2*Ro2*^Rol )/2O0 
Al = (Ro2-Rol-900*A2 )/10 
A0=Rol-40*Al-1600*A2 
Ro*A0*^Egr* ( A1 +Egr*A2 ) 

RETURN Ro 
FNEND 

DEF FNPreg( T , Egr ) 

Pr = FNCpeg( T,Egr )*FNNueg< T , Egr )*FNRoeg( T,Egr )/FNKeg(T,Egr ) 

RETURN Pr 
FNEND 

DEF FNKeg< T c , Egr ) 

RANGE OF VALIDITY: -20 TO 20 DEG C 

Tk*Tc+273.15 

K 1 =2. 2824708E- 1 +T^ *( 5. 5989286E-4 + Tk *3.57142866-7 ) 

K2-2. 58466 16E-1+Tk*( 2. 397857 lE-4 + Tk *7. 142857 lE-7 > 

K3 = 3.2138932E-1-Tk*< 3.0042857E-4-Tk*l .4285714E-6 ) 

A2'(K3-2*K2 + K1 )/200 
A1 = (K2-K1-900*A2 )/10 
A0=K1-40*A1-1600*A2 
K-A0 + Egr* < A1 +Egr*A2 ) 
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2E1C 


RETURN K 




:bi5 


FNEND 




:b:o 


OEF FNTanh( / ) 




2B25 


P = EXP( X ) 




2630 


Q*l/P 




2B55 


Tanh=^P-Q)/(P+Q> 




2640 


RETURN Tann 




2645 


FNEND 




2650 


DEF FNTvsvkU ) 




2655 


COM /Cc/ C(7) 




2660 


T*C<0) 




2665 


FOR 1=1 TO 7 




2670 


T=T*^C( I 




2675 


NEXT I 




26B0 


RETURN T 




2685 


FNEND 




2690 


DEF FNBeta(T ) 




2695 


Rop»FNRho(T+.l ) 




2700 


Ron=FNRho< T- . 1 ) 




2705 


Beta=-2/<Rop+Ron)*(Rop-Rom)/ 


t n 


2710 


RETURN Beta 




2715 


FNEND 




2720 


DEF FNPsaUTc ) 




2725 1 


0 TO 80 deg F CURVE FIT OF Psat 


2730 


Tf=l ,B«Tc+32 




2735 


Pa=5.945525*Tf«( . 1 5352082 + T f • ( 1 .4840963E-3+Tf *9.81506718-8 ) ) 


2740 


Pg*Pa- 1 4 . 7 




2745 


IF Pg>0 THEN « +-P5IG,- 


-“in Hg 


2750 


Psat=Pg 




2755 


ELSE 




2760 


P5at=Pg*29.92/14.7 




2765 


END IF 




2770 


RETURN Psat 




( /b 


r ixc-ivu 




27B0 


DEF FNFncaK I ,T,Fr) ) 




2785 


DIM Bl<4 ),B2(4),B3(4 ),M1(4) 


,M2<4 ),M3<4 ) 


2790 


DATA 3.48835E-3,6 .71749E-5, 


-5.8103E-3,-5.5079E-3.4. 125E-4 


2795 


DATA -1 .47B21E-2,2.53207E-3 


,-2.7539BE-3,-4.309l3E-3,-3.10937£- 


2796 


DATA 1.04BE-3,1 .01292E-2,1. 


114E-2,-3.86717E-3,0.0 


2B00 


DATA 1 .71B25E-3,2.70428E-3, 


1 .93749E-3,! .92237E-3, 2.566 4BE-3 


2805 


DATA 2.034B2E-3,2.B6716E-3, 


2. 095 17E-3, 2. 124336-3,2. 93127E-5 


2B0B 


DATA 1 .7B81E-3,2.70622E-3,2 


.00121E-3, 2. 20406-3,0.0 


2810 


READ BK * ),B2< ♦ ),B3( * ),M1( ♦ 


),M2( ♦ ),M3( * ) 


2811 


IF T<-7 THEN 




2815 


Zl-BK I ) + Ml( I )*Fn 




2B20 


Z2 = B2< I ) + M2< I )*Fn 




2B21 


5 = < Z2-Z1 )/13 




2822 


C=Zl+5*20 




2824 


ELSE 




2825 


Z1*B2< I )+M2< I )*Fn 




2B2B 


Z2 = B3< I )+M3< I )*Fn 




2B27 


S=( Z2-Z1 )/16 




2828 


C=Z1+S*7 




2830 


END IF 




2B35 


Mdot=C+S*T 




2B40 


RETURN Mdot 




2B45 


FNEND 




2850 


SUB Xyread 




2B55 


BEEP 
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2960 

2665 

2570 

2676 

2880 

2885 

2890 

2595 

2900 

2905 

2910 

2915 

2920 

2925 

2930 

2935 

2940 

2945 

2950 

2955 

2990 

2965 

2970 

2975 

2980 

2985 

2990 

2995 

3000 

3005 

3010 

3015 

3020 

3025 

3030 

3035 

3040 

3045 

3050 

3055 

3065 

3070 

3075 

3080 

3085 

3090 

3095 

3100 

3105 

3110 

3115 

3120 

3125 

3130 

3135 

3140 

3145 

3150 

3155 



INPUT “ENTER FILE NAME’‘,File$ 

BEEF 

INPUT "ENTER NUMBER OF X,Y PAIRS", N 
ASSIGN ©File TO FileS 
FOR 1=1 TO N 

ENTER @File.X,Y 
PRINT X,Y 
NEXT I 
SUBEND 
SUB Purge 
BEEP 

INPUT "ENTER FILE NAME TO BE DELETED ", F 1 1 e$ 

PURGE Files 
GOTO 2910 
SUBEND 
SUB Wilson 

COM /Will/ Doa(4),Dia(4),Kna(4),Iact 

COM /Wil2/ Delta, Isat ,N5ets,Hod,Cia( 3 ) , Alpaa( 3 ) 

DIM Fn(4 ),Enf(8 ),Tp(3),T(8 ),Xa(20),Ya(20 ) 

BEEP 

INPUT "PLEASE RE-ENTER NAME OF FILE",File$ 

ASSIGN ©File TO FileS 

INPUT "ENTER TUBE CODE",Itube 

BEEF 

INPUT “GIUE A NAME FOR XY FILE",Xy$ 

CREATE BOAT Xy$,5 
ASSIGN ©Xy TO XyS 
L=1 .2192 
Do=Doa( I tube ) 

Di®Dia( I tube ) 

Kn*Kna ( 1 1 ube ) 

Ax=PI *Di "2/4 I Cross-sectional area 
Ao=PI *Do*L 

Rn“Do*LOG( Do/Di )/(2*Kn) 

Initial values 
Tf=Tsat 
Alpa=.E55 
Ci=5. 172 
6=9.81 

i Deg=fe3 

Iend=5 ICHANGE TO 4, IF FIUE TUBES IN BUNDLE 
IF Iact<10 THEN 
Ibeg=Iac t 
Iend= I ac t 
END IF 

FOR I»Ibeg TO lend 
Sx = 0 
Sy=0 
Sxs = 0 
Sxy=0 
Jset =0 

ASSIGN ©File TO FileS 

ENTER ©File;DtgS,Itube,Egrat,Ddl,Dd2,Dd3,Dd4,Dd5 
ENTER ©Fi 1 e i Fn( ♦ ) ,Enf ( * ) , Tp( * ) 

FOR J=0 TO 8 

T( J )*FNTvsv(Enf ( J ) ) 

NEXT J 

Tvap=( T( 0 ) + TC 1 ) + T( 2 ) )/3 
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21E0 

31B5 

3170 

3175 

3180 

3185 

3190 

3195 

3300 

3205 ‘ 

3210 ' 

3215 

5220 

3230 

3235 

3240 

3245 * 

3250 

3255 

32E0 

32E5 

3270 

3275 

3280 

3285 

3290 

3295 ‘ 

3300 I 

3305 

3510 

3315 

3320 

3325 ‘ 

3330 

3335 

3340 

5345 

3350 

3355 

3350 I 

33E5 

3370 

3375 

3380 

3385 

3395 

3400 

5405 

3410 

3415 

3420 

3425 

3430 

3435 

3440 

3445 

3450 

3455 I 



T1:q*( T f 3 4 ) )/2 

IF lsat*0 THEN 
I s a t = T 1 1 q 
ELSE 

T sat = Tvap 
END IF 

Grad = FNGrao( T( I +5 ) ) 

Delt=ABS(Tp( I )/(Grad*10) ) 

Tavg*T( 1 + 5 ) + Delt*.5 

Water/Et hy lane Glycol Mixture Properties 
Rhoeg*FNRhoeg( Tavg , Egrat ) 

Nueg=FNNueg( Tavg , Egrax ) 

Mueg*Nueg*Rhoeg 
Cpeg=FNCpeg( Ta vg , Egra t ) 

Keg=FNKeg ( Ta vg , Egr a t ) 

Preg*Cpeg*Mueg/Keg 

Mdot=FNFmcal( I,T( 1+5 ),Fn( I ) ) 

Ueg^Mdot /( Rhoeg* Ax ) 

Reeg-Veg*Di /Nueg 

Pes*4*Mdot/(PI* Mu eg*(Di-4* Delta ) ) 

Odot=Mdot+Cpeg*Delt 

Odp-Odot /Ao 

Lntd«Delt/L0G( ( Tsat-T( 1+5 ) )/( Tsat-T(I + 5 )-Delt ) ) 
Uo*0dp/Lmt d 

0nega = ( 1 +5 . 484E-3*Preg ' . 7+ ( Res/Hod ) " 1 . 25 )" . 5 

R-114 Properties 
Hfg«FNHf g( Tsat ) 

Kf*FNK(Tf ) 

Rhof*FNRho( Tf ) 

Muf-FNMu( Tf ) 

F«( Kf "3 + Rhof ''2*G + Hfg/(Muf *Do*0dp ) )" .33333 

Ho*Alpa*F 

Two*Tsat “Odp/Ho 

Tf“Tsat/3+2*Two/3 

Y= ( 1 /Uo-Rm )*F 

X*Do*F/ ( Keg ♦Onega ) 

PRINT "0MEGA=" ; Onega ; "F=" ;F; "X*" ;X; ”Y=" : Y 
Xa(Jset)=X I INEFFICIENT (MODIFY LATER) 

Ya( Jset )*Y 
Sx=Sx+X 
5y*5y+Y 
Sx s = Sx s + X *X 
J set - Jset + 1 

IF Jset^Nsets THEN 3135 
ASSIGN @Fiie TO ♦ 

Slope*(N3ets*Sxy“Sx*Sy )/( Nsets*Sxs-Sx''2 ) 

Intcpt=(Sy-Slope*Sx )/Nsets 

Cl c«l /Slope 

Alpac*l / Int cpt 

Cerr=ABS( (Ci-Cic )/Cic ) 

Aerr*ABS( (Alpac-Alpa )/Alpac ) 

IF Cerr>.001 OR Aerr>.001 THEN 
Alpa = (Alpa + Alpac )*.5 
Cl “( C 1 +Ci c )♦ . 5 

PRINT ''CIC=’‘ ; Cl c ; ’‘ALPA-“ ; Alpa 
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54BC 

346E 

3470 

3475 

54B0 

3485 

3490 

3495 

3500 

5505 

3510 

3515 

3520 

3525 

3530 

3535 

3540 

3545 

3550 

3555 

3560 

5565 

35"’0 

3575 

3580 

3565 

3590 

3555 

3600 

3605 

3610 

3615 

3620 

3625 

3630 

3655 

3640 

3645 

3650 

3655 

3660 

3665 

3670 

3675 

3680 

3685 

3690 

3695 

3700 

3705 

5710 

3715 

3725 

3730 

3735 

3740 

3745 

3750 

3755 



GOTO 3100 

END 

BEEP 

6EEF 

PPINTEF 15 1 

PRINT "CIC=" ;Cic; .Pilpa 

Ci a ( I >=C 1 c 

Alpaa^ I )=f‘tlpac 

PRINTER IS 701 

FOR J*0 TO Nsets-l 

OUTPUT @Xy. Aa( J ),Ya( J ) 

NEXT J 

PRINTER IS 1 
NEXT I 

ASSIGN @Xy TO * 

SUBEND 
SUB Nusselt 

CON /Nus/ I in, Tsat , Qdp , Hoc 
Do=.0155 
ho- 1 000 
IF 1 1 n-0 Then 
BEEF 

INPUT "ENTER TSAT AND HEAT FLUX ", T = at , Qdp 
END IF 

Hf g = FNHf g( Tsat ) 

Two=T5at-Qdp /Ho 
Tf=Tsat/3+2*Two/3 
Kf*FNK(Tf ) 

Rho f =FNRho( T f ) 

Nuf=FNNu( Tf > 

Hoc=.655*( Kf "3*Rhof ".2*9.81*Hfg/(Muf*0c*0dp ) )". 335333 
IF ABS((Ho-Hoc )/Hoc )>.001 THEN 
Ho® ( Ho+Hoc )♦ . 5 
GOTO 3585 
END IF 

IF lin = 0 THEN PRINT "HO®" '.Hoc 

SUBEND 

SUB Flotl 

DIM >aa(4 ) 

PRINTER IS 705 

Idv=l 

BEEP 

input "OK TO USE DEFAULT UALUE5 ( 1 =DEF = Y , 0®N ) " , I dv 
IF Idv-1 THEN 
Itn=2 
Xni n= 1 
Xnax®5 
Xstep=l 
Ymn=0 
Ynax«2 . 0 
Ys tep® . 5 
ELSE 

INPUT "ENTER MINIMUM AND MAXIMUM X-UALUES " , Xrnin Xnar. 
BEEP 

INPUT "ENTER MINIMUM AND MAXIMUM Y-UALUES " , Yt^in 'ins>. 
BEEP 

INPUT "ENTER STEP SIZE FOR X-AX I S " , X s t ep 
BEEP 

INPUT "ENTER STEP SIZE FOR Y-AX IS" , Ystep 
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3750 

Z'"65 

37^0 

3775 

3760 

5785 

3790 

57S5 

5800 

5805 

5810 

3815 

5820 

5525 

3830 

3835 

3840 

3845 

3650 

3855 

3850 

3885 

3570 

3875 

3880 

3885 

3890 

3895 

3900 

3905 

3910 

3915 

3920 

3925 

3930 

3935 

3940 

3945 

3950 

3955 

3360 

3965 

3970 

3975 

3980 

3985 

3990 

3995 

4000 

4005 

4010 

4015 

4020 

4025 

4030 

4035 

4040 

4045 



END IF 

PRINT ‘'IN.SPlilP 2300, 1800,8300,6800. 

PRINT "SC 0, 100,0, 100. TL 2.0." 

5fx = 100/( XnaA-Xnin ) 

Sfy = 100/(Yna>-Yrun) 

PRINT "PU 0,0 PD" 

FOR Xa=Xnin TO Xnax STEP Xstep 
X = ( Xa-Xnin )*5f / 

PRINT "PA";X,",0. XT. " 

NEXT Xa 

PRINT "PA 100,0. PU; " 

PRINT "PU PA 0,0 PD" 

FOR Ya=Vnin TO Ynax STEF Ystep 

Y = (Ya-Ynin )*5fy 
PRINT "PA 0, ";Y. "YT" 

NEXT Ya 

PRINT "PA 0, 100 TL 0 2" 

FOR Xa=Xmin TO Xnax STEP Xstep 
X = ( Xa-Xnin ;*Sf x 
PRINT "PA" . X, " , 100. XT- 
NEXT Xa 

PRINT "PA 100,100 PU PA 100,0 PC" 

FOR Ya=Ymin TO Ymax STEP Ystep 

Y = ( Ya->nin )*Sf y 

PRINT "PD PA 100, " ,Y, "YT" 

NEXT Ya 

PRINT "PA 100,100 PU" 

PRINT "PA 0,-2 SR 1.5,2" 

FOR Xa*Xnin TO Xnax STEP Xstep 
X==( Xa-Xn : n >*S f x 
PRINT "PA" ;X, " ,0; " 

PRINT “CP -2,-l;LB";Xa: "" 

NEXT Xa 

PRINT "PU PA 0,0" 

FOR Ye=Ynin TO Ynax STEP Ystep 

Y = ( Ya-Yn 1 n >*S f y 
PRINT "PA 0, " ; Y, "" 

PRINT "CP -4,-.25:LB" :Ya; "" 

NEXT Ya 
BEEP 

INPUT "SELECT MODE ( 0=HN/H1 , 1 =HN( a vg ) /H 1 > ",I sr^ 

I sn = I sn+ 1 
IF Idv=l THEN 

IF IsM=l THEN Ylabel$="HN/Hl " 

IF Isn=2 THEN Ylabel$="HN( avg )/Hl" 
Xiabel$="Tube Nunber" 

ELSE 

BEEP 

INPUT "ENTER X-LABEL" , XlabeiS 
BEEP 

INPUT "ENTER Y-LABEL" , YlabelS 
END IF 

PRINT “SR 1.5,2;PU PA 50,-10 CP -LEN( X labels )/2 ; 
PRINT "PA -11,50 CP 0,“;-LEN<Ylabel$>/2*5/6; "01 0 
PRINT "CP 0,0 DI" 

Okp = l 
BEEP 

INPUT "WANT TO PLOT DATA FROM A FILE (1=DEF*Y,0=N 



"0;LB" iXlabelS; "" 
, 1 ;LB" ;Ylabel$; " " 



)7" ,Okp 
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4055 

4OB0 

40E5 

4070 

4075 

4080 

4085 

4090 

4095 

4100 

4105 

4110 

4115 

4120 

4125 

4130 

4135 

4140 

4145 

4150 

4155 

41B0 

4155 

4170 

4175 

4180 

4185 

4190 

4195 

4200 

4205 

4210 

4215 

4220 

4225 

4230 

4235 

4240 

4245 

4250 

4255 

4250 

4255 

4270 

4275 

4280 

4285 

4290 

4295 

4300 

4305 

4310 

4315 

4320 

4325 

4330 

4335 

4340 

4345 

4350 



BEEF 

INPUT •'ENTER THE NAME OF THE DATA FILE",DfileS 

ASSIGN @Filft TO DfileS 

BEEP 

INPUT -ENTER THE BEGINNING RUN NUMBER”, Md 
BEEF 

INPUT ‘'ENTER THE NUMBER OF X-V PAIRS STORED ", Nset s 
BEEP 

INPUT "SELECT A SYMBOL FOR THE PLOTTER <1 = * , 2 = + , 3=c . 4 = o , 5* ^ , 5 > 

PRINT "PU DI" 

IF 5y=l THEN PRINT "SM*" 

IF 5y*2 THEN PRINT "SM+" 

IF Sy=3 THEN PRINT "SMc" 

IF 5y=4 THEN PRINT "SMo" 

IF 5y*5 THEN PRINT "SM''" 

FOR 1*1 TO Nsets 
FOR J=0 TO 3 

ENTER @File: Vaa( J ),D 
IF J*0 THEN Ytop=>aa(0) 

Yaa( J )*Vaa( J )/Ytop 
NEXT J 

FOR J*0 TO 3 

X*( J + l-Xmin )*Sf p- 
Y = ( Yaa( J )-Ynin )*Sf y 
PRINT "PA",X,Y,"" 

NEXT J 
NEXT I 
BEEP 

ASSIGN @File TO * 

GOTO 4040 
END IF 
BEEP 

INPUT "LIKE TO PLOT THE NU5SELT RELATION ( 1*Y, 0*N )•? " , Oknus 
PRINT "PU:5M" 

IF Oknus*l THEN 

FOP Xa*Xnin TO Xna^ STEP Xstep/50 
X=(Xa-Xnin)*SfA 

IF Isn=l AND Xa^'Xnin THEN Ya = Xa . 75-( Xa-1 ) . 75 
IF Isn=2 AND Xa>Xnin THEN Ya*Xa'‘ ( - . 25 ) 

IF Xa=Xnin THEN Ya*l 
Y=( Ya-Ynin )*Sf y 
PRINT "PA" ,X, Y, "PD" 

NEXT Xa 
BEEP 

PRINT "PU’‘ 

INPUT "MOVE THE PEN TO LABEL THE NU5SELT LINE", Ok 
PRINT "LBNusselt" 

END IF 

IF Isn*2 THEN 
BEEP 

INPUT "LIKE TO PLOT EXPTL CURUE ( 1 =Y , 0*N ) " , Okex 
Nq*0 

IF Okex*l THEN 
BEEP 

INPUT ‘’ENTER THE EXPONENT", Ex 
FOR Xa*Xnin TO Xnax STEP Xstep/10 
Nq»Nq+l 
Ya-Xa"(“Ex ) 

X*( Xa-Xmin )*Sf x 
Y*( Ya-Ynin )*Sf y 
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4355 

43B0 

43B5 

437C 

4375 

4 

4385 

4390 

4395 

4400 

4405 

4410 

4415 

4420 

4425 

4430 

4435 

4440 

4445 

4450 

4455 

4 460 

44B5 

4470 

4475 

4480 

4485 

4490 

4495 

4500 

4505 

4510 

4515 

4520 

4525 

4530 

4535 

4540 

4545 

4550 

4555 

45G0 

4565 

4570 

4575 

4580 

4585 

4590 

4595 

4B00 

4505 

4510 

4B15 

4520 

4525 

4550 

4535 

4540 

4545 



IF Nq MOD 2=0 THEN 
PRINT "PA’‘,X,Y/‘PD" 

ELSE 

PRINT ”PA" ,X ,Y, ‘’PU“ 

END IF 

NtA I A a 
PRINT "PU" 

BEEP 

INPUT "MOUE PEN TO LABEL AND HIT ENTER", 0^ 

PRINT "LBs=0" 

PRINT "PR -1 0" 

PRINT "LB":Ex;"" 

GOTO 4300 
END IF 
END IF 
GOTO 4510 
BEEP 

INPUT "LIKE TO PLOT KERN RELATIONSHIP (1 =Y , 0=N )7 " , Yes 
IF Yes=l THEN 

FOR Xd=Xnin TO Xnax STEP Xstep/20 
Ya = xa'*( -1/5 ) 

X = ( Xa-Xnin )*Sf . 

Y = ( Ya-Vnin )*Sfy 
PRINT “PA" ,X , Y/'PD" 

NEXT Xa 
PRINT "PU" 

BEEP 

INPUT "MOUE THE PEN TO LABEL KERN RELATI ONSHIP , Ok 
PRINT "LBKern;PU" 

END IF 

PRINT "PU PA 0,0" 

BEEP 

INPUT "LIKE TO PLOT EISSENBERG RELATION (1 =Y , 0=N )7 " , Oke i 
IF Okei=l THEN 

FOR Xa=Xnin TO XnaN STEP Xstep/10 
Ya=.E^.42*Xa^(-.25) 

X = ( Xa-Xnin )«Sf>. 

Y=( Ya-Ynin )*Sfy 
PRINT "PA" ,X ,Y, "PD” 

NEXT Xa 
PRINT "FU" 

BEEF 

INPUT "MOUE THE PEN TO LABEL THE EISSENBERG LINE", Ok 
PRINT "LBEi ssenberg ; PU " 

END IF 

PRINT "PU SP0" 

SUBEND 

DEF FNPvst (Tc ) 

COM /Fid/ Ift 
DIM K( 8 ) 

IF lft*0 THEN 
BEEP 

PRINT "PUST CORRELATION NOT AVAILABLE FOR R-114" 

STOP 
END IF 

IF Ift*l THEN 

DATA -7.591234554,-25 .0802359B,-1B8. 170B54B, 54. 23285504, -11 8. 9545225 
DATA 4.15711732,20.9750575,1.59,5 
READ K( ♦ ) 
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4B50 

4B55 

4BB0 

4BB5 

4B70 

4B75 

4B80 

4B85 

4B90 

4B95 

4700 

4705 

4715 

4720 

4725 

4730 

4735 

4740 

4745 

4750 

4755 

4750 

4755 

4770 

4775 

4780 

4785 

4790 

4795 

4800 

4805 

4810 

4815 

4820 

4825 

4830 

4835 

4840 

4845 

4850 

4855 

4850 

4855 

4870 

4875 

4880 

4885 

4890 

4895 

4900 

4905 

4910 

4915 

4920 

4925 

4930 

4935 

4940 

4945 



T=(Tc+273.15)/647.3 

5urv=0 

FOR N*0 TO 4 

Sum*Sufn+K(N)*( 1-T >''(N+1 ) 

NEXT N 

Br-5un/( T*( 1+K( 5 )* (1-T ) + K( 5 )♦( 1 -T >"2 ) )-(l-T )/( K ( 7 )♦ (1-T )"2 + K( 8 ) ) 
Pr-EXP(Br) 

P«22l20000*Pr 
END IF 

IF Ift*2 THEN 

Tf*Tc*l .8+32+459.5 

P*=10"( 33. 0555-4350. 98/Tf -9. 2635*LGT(Tf ) + 2 . 0539E-3*T f > 

END IF 

IF Ift*3 THEN 

A-9.394B85-3055.1/(Tc + 273. 15) 

P*133.32*10''A 
END IF 
RETURN P 
FNEND 

DEF FNHfg(T ) 

COM /Fid/ Ift 
IF Ift-0 THEN 
Tf=T*l .8+32 

Hfg-5. 1451 558E+1-Tf*( 5.951 079E-2 + Tf*( 1 . 3988588E-4+1 . 9507843E-7*T f ) ) 
Hfg=HfQ*2325 
END IF 

IF Ift=l THEN 

Hfg*2477200-2450*(T-10 ) 

END IF 

IF Ift-2 THEN 
Tf*T*1.8+32 

Hfg=7.0557857E+l-Tf*( 4.838052E-2+1 .2619048E-4*Tf ) 

Hfg*Hfg*2325 
END IF 

IF Ift-3 THEN 
Tk-T+273.15 

Hfg=l .35254E+B-Tk*( 5. 38253E+2+TM .747452 ) 

END IF 
RETURN Hfg 
FNEND 

DEF FNMu( T ) 

COM /Fid/ Ift 
IF lft=0 THEN 
Tk»T+273.15 

Mu«EXP( -4.4535+1011 .47/Tk )M .E-3 
END IF 

IF Ift=l THEN 

A*247.8/(T+133.15) 

Mu*2.4E-5M0"A 
END IF 

IF Ift«2 THEN 

riu-8.9529819E-4-T*( 1 . 1094509E-5-T*5 . 565 829E-8 > 

END IF 

IF Ift=3 THEN 

Tk»l/(T+273.15) 

Mu»EXP( -11.0179+Tk*( 1 .744E+3-Tk*( 2.80335E+5-TkM.12551E+8> ) ) 

END IF 
RETURN Mu 
FNEND 
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4950 

4955 

4960 

4965 

4970 

4975 

4980 

4985 

4990 

4995 

5000 

5005 

5010 

5015 

5020 

5025 

5030 

5035 

b04^ 

5045 

5050 

5055 

5060 

5065 

5070 

5075 

5080 

5085 

5090 

5095 

5100 

5105 

5110 

5115 

5120 

5125 

5130 

5135 

5140 

5145 

5150 

5155 

5160 

5165 

5170 

5175 

5180 

5185 

5190 

5195 

5200 

5205 

5210 

5215 

5220 

5225 

5230 

5235 

5240 



DEF FNUvst( It ) 

COM /Fid/ Ift 
IF lft*0 THEN 
8EEP 

PRINT “UUST CORRELATION NOT AC'AILAELE FOR R-114" 

STOP 
END IF 

IF Ift»l THEN 
P-FNPvst (Tt ) 

T-Tt+273.15 

X»1500/T 

Fl*l/( 1+T*1 .E-4 ) 

F2«( l-EXP(-X) )"2.5*EXP(X >/X".5 
0015*F1 000942 *F2-. 0004882 *X 
K-2*P/< 461 .52*T ) 

l + ( 1+2*B*K )" . 5 )/K 

END IF 

IF Ift=2 THEN 
I f » I t ♦ 1 . 

l>=13. 955357-Tf*( . 16127262-Tf *5 . 1 7261 90E-4 ) 

»;*U/16.018 
END IF 

IF Ift*3 THEN 
Tk*Tt+273. 15 
P-FNPvst ( Tt ) 

U*133.95*Tk/P 
END IF 
RETURN U 
FNEND 

DEF FNCp(T) 

COM /Fid/ Ift 
IF lft=0 THEN 
Tk»T+273. 15 

Cp*.40118+Tk*( 1 .65007E-3^Tk*( 1 . 51494E-6-Tk *6 . 67853E- 10 ) ) 

END IF 

IF Ift*l THEN 

Cp*4.21120858-T*(2.2682GE-3-T*(4.42361E-5+2.71428E-7*T ) ) 

END IF 

IF Ift=2 THEN 

Cp*9.2507273£-1+T*( 9.3400433E-4+1 .7207792£-6*T ) 

END IF 

IF Ift-3 THEN 
Tk*T^273.15 

Cp*4. 1868*( 1 .6884E-2 + TkM 3.35083E-3-Tk*(7.224E-6-Tk*7.6l748E-9 ) ) ) 
END IF 

RETURN Cp*1000 
FNEND 

DEF FNRho(T) 

COM /Fid/ Ift 
IF lft*0 THEN 
Tk*T-^273.15 
X=l-(1.8»Tk/753.95) 

Ro*36. 32 + 61 . 1464 14*X"( 1/3 ) + 16 . 4 1 8015*X+ 17 . 476838»X '' .5+1 . 119828 *X" 2 
Ro=Ro/ . 062428 
END IF 

IF Ift=l THEN 

RO-999.52946+TM . 01269-T* ( 5 . 482513E-3-T* 1 . 234147E-5 ) ) 

END IF 

IF Ift-2 THEN 
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5Z45 

5Z50 

5255 

52E0 

52G5 

5270 

5275 

5280 

5285 

5290 

5295 

5300 

5305 

5310 

5315 

5320 

5325 

5330 

5335 

5340 

5345 

5350 

5355 

5380 

5385 

5375 

5580 

5385 

5390 

5395 

5400 

5405 

5410 

5415 

5420 

5425 

5430 

5435 

5440 

5445 

5450 

5455 

5480 

5485 

5470 

5475 

5480 

5485 

5490 

5495 

5500 

5505 

5510 

5515 

5520 

5525 

5530 

5535 

5540 



Ro=l .G20747 9E + 3-T*( 2 .2185348 + 1*2.3578291 £-3 ) 

END IF 

IF Ift=3 THEN 

Tk=T+273. 15-338. 15 

Uf»9.24848E-4 + TI' * ( 8 . 2796E-7 + Tk * ( 9 . 2444E- 1 0 + TK *3 . 057E-1 2 ) ) 
Ro-l/'Jf 
END IF 
RETURN Ro 
FNEND 

DEF FNPr( T > 

Pr = FNCp{ T )*FNMu( T )/FNK(T ) 

RETURN Pr 
FNEND 

DEF FNK(T ) 

COM /Fid/ Ift 

IF lft“0 THEN K'.071-.0002G1*T 
IF Ift=l THEN 

X«(T + 273. 15 )/273. 15 

K = -.92247 + X*( 2.8395-X*{ 1 .8007-X*( . 52577- . 07344 *X ) ) ) 

END IF 

IF Ift=2 THEN 

K=8.2095238E-2-T*( 2.221 428GE-4+T*2.3809524E-8 ) 

END IF 

IF Ift*3 THEN 
Tk-T+273. 15 
END IF 
RETURN K 
FNEND 

DEF FNHf(T) 

COM /Fid/ Ift 
IF Ift-0 THEN 
BEEP 

PRINT “HF CORRELATION NOT FOR R-114" 

STOP 
END IF 

IF Ift=l THEN 

Hf«T*( 4 .203849-T* ( B. 88132E-4-T*4 .551G0317E-6 ) ) 

END IF 

IF Ift=2 THEN 
Tf=T*l .8+32 

H f = 8. 207857 1+Tf*( . 1 94G7857 + T f * 1 .321428GE-4 ) 

Hf=Hf*2.32G 
END IF 

IF Ift=3 THEN 

Hf=250 t TO BE UERIFIED 
END IF 

RETURN Hf*1000 

FNEND 

SUB Plot2 

COM /Drl/ S t ar , Syn , Icon 

COM /Fid/ Ift 

DIM C(9),Xya(7),Doa(3 ) 

DATA 0.0158,0.0158,0.0158,0.0158 
READ Doa( * ) 

Fui*=l 

PRINTER IS 1 
BEEP 

PRINT USING "4X "Select Option X-Y Limits **"" 

PRINT USING "GX,“‘'0 Use default values""" 
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5545 

5550 

5555 

5550 

5565 

5570 

5575 

5580 

5565 

5590 

5595 

5600 

5605 

5610 

5615 

5620 

5625 

5630 

5635 

5640 

5645 

5650 

5655 

5660 

5665 

5670 

5675 

5680 

5685 

5690 

5695 

b /0W 

5705 

5710 

5715 

5720 

5725 

5730 

5735 

5740 

5745 

5750 

5755 

5760 

5765 

5770 

5775 

5780 

5785 

5790 

5795 

5800 

5805 

5810 

5815 

5820 

5825 

5650 

5835 

5840 



FPIN"^ USING Use neui values 

INPUT Cr d 
BEEF 

INPUT "ENTEF TUBE CODE'Mcode 
Do*Dca( 1 1 ube ) 

Iht=2 

BEEP 

PRINT USING "4X , ""Select option: 

PRINT USING "6X,""0 h versus q 

PRINT USING "6X/’*‘l q versus Delta-T 

PRINT USING "6X,’’"2 h versus Deita-T (default) 

INPUT Iht 
PRINTER 15 705 

IF OKd=0 THEN lAXIS DEFAULT UALUES 
IF lht»0 THEN < ( h vs q ) 

Ym n=*0 

Y n a X = 6 0 
Ys t ep = 1 0 
Xmin* . 2 
Xnax^l .4 
Xstep- . 2 

END IF 

IF Iht-1 THEN !<q vs t ) 

Xni n=0 
Ym n*0 
Ynax= . 5 
Xnax * 1 5 
X s t ep®3 
Ystep= . 1 
END IF 

IF Iht»2 THEN ' < h vs t ) 

Xn 1 n*U 
Yni n=0 
Xnax=50 
Ynax =6 
Xs t ep- 1 0 

Y s t e p = I 
END IF 

END IP 

IF Okd=l THEN 
BEEP 

INPUT "ENTER MINIMUM AND MAXIMUM X-UALUES " , Xni n , Xnax 
BEEP 

INPUT "ENTER MINIMUM AND MAXIMUM Y-UALUES " , Yni n , Ynax 
BEEP 

INPUT "ENTER STEF SIZE FOR X-AX I S " , X s t ep 
BEEP 

INPUT "ENTER STEP SIZE FOR Y-AX IS " , Ystep 
END IF 
BEEF 

PRINT "IN;SP1;IP 2300,1800,8300.6800;" 

PRINT "SC 0, 100,0, 100;TL 2,0;" 

Sfx=100/( Xnax-Xnin ) 

Sfy-100/(Ynax-Ynin ) 

BEEP 
I CQ = 0 

INPUT "LIKE TO BY-PASS CAGE ( 1 =Y , 0-N»DEFAULT ) 7 " , I cq 
IF IcQ-1 THEN 6175 
PRINT "PU 0,0 PD" 

FOR Xa=Xnin TO Xna^ STEF Xstep 
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5845 X = ( Xa-Xnin )»S f > 

5850 PRINT "Pf^” : X , ■ ,0. XT : 

5855 NEXT Xa 

5660 PRINT -FiS 1O0,C,PU," 

5665 PRINT "PU PA 0,0 PD" 

5570 FOR Ya = Ynin TO Yna;-. STEP Ystep 
5875 Y = ( Ya-Yrnin )‘Sf y 

5660 PRINT "PA 0, " . Y, "YT” 

5885 NEXT Ya 

5890 PRINT "PA 0,100 TL 0 2" 

5895 FOR Xa*Xnin TO Xmax STEP Xstep 

5900 X=( Xa-Xnin >*Sfx 

5905 PRINT "PA";X,",100; XT" 

5910 NEXT Xa 

5915 PRINT "PA 100,100 PU PA 100,0 PD" 

5920 FOR Ya=Yrun TO Ynax STEP Ystep 

5925 Y = ( Ya-Ynin )»Sfy 

5930 PRINT "PD PA 100,",Y,"YT" 

5935 NEXT Ya 

5940 PRINT "PA 100,100 PU" 

5945 PRINT "PA 0,-2 SR 1.5,2" 

5950 FOR Xa=Xnin TO Xnax STEP Xstep 

5955 X=(Xa-Xnin)*Sfx 

5960 PRINT "PA";X,",0;" 

5965 IF Xa<l AND Xa<>0 THEN PRINT "CP -1 . 5 , - 1 ; LB0 i PR -1 , 0 : LB “ . Xa " 

5970 IF Xa=0 THEN PRINT “CP -.5,-1 :LB0" 

5975 Xin=0 

5980 IF Xa MOD 1-0 THEN Xin-1 

5985 IF Xa>-10 THEN PRINT "CP -2 , -1 ; LB " ; Xa ; " " 

5990 IF Xa>l AND Xa<10 AND Xin-1 THEN PRINT "CP -1 . 25 , -1 ; LB " ; Xa i " " 

5995 IF Xa>l AND Xin-0.THEN PRINT "CP -2 , -1 ; LB" : Xa ; " " 

6000 IF Xa-1 THEN PRINT "CP -1,-1;LB1.0" 

6005 NEXT Xa 

6010 PRINT "PU FA 0,0" 

6015 FOR Ya-Ynin TO Ynax STEP Ystep 

6020 Y=( Ya-Yrun )*Sfy 

6025 PRINT "PA 0, " . Y. " " 

6035 IF Ya<l AND Ya^>0 THEN PRINT "CP -4 , - . 25 ; LB0 ; PR -2 , 0 . LB " ; Ya ; ” " 

6040 IF Ya-0 THEN PRINT "CP -2,-.25.LB0" 

6045 IF Ya>l AND Iht<2 THEN PRINT "CP -5 , - . 25 . LB " ; Ya : " " 

6050 IF Ya>9 AND Iht-2 THEN PRINT "CP -4 , - . 25 ^ LE " ; Ya ; " " 

6055 IF Ya>0 AND Ya<10 AND Iht-2 THEN PRINT "CP -3 . - . 25 , L8" ; Ya ; " " 

6060 I IF Ya-1 THEN PRINT "CP -4 , - . 25 ; LB 1 . 0 " 

6065 NEXT Ya 

6070 IF Okd=l THEN 

6075 BEEP 

6080 INPUT "ENTER X-LABEL" , XlabeU 

6085 BEEP 

6090 INPUT "ENTER Y-LABEL " , Y1 abe 1 $ 

6095 END IF 

6100 IF IhtOl THEN 

6105 PRINT "SR 1.5,2iPU PA -12,35;DI 0,l;LBh;PR l,0.5;LBoiPR -l,0.5iLB/( 

kU/n" 

6110 PRINT "PR -1,0.5;SR 1.1.5;LB2;SR 1.5,2:PR .5,.5;LB.iPR .5,0iLBK>" 

6115 ELSE 

6120 PRINT "PA -12,39iDI 0 , 1 ; LBq/ ( MW/n ; SR 1,1.5;PR - 1 , 0 . 5 . LB2 i SR 1.5,2.P 

R 1.0, 0.5; LB)” 

6125 END IF 

6130 IF Iht-0 THEN 
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6135 

1.5,: 

6140 

6145 

:,c 

6150 

6155 

6160 

61E5 

6170 

6175 

6180 

6185 

6190 

6195 

6200 

6205 

6210 

6215 

6220 

6225 

6230 

6235 

6240 

6245 

6250 

6255 

6260 

6265 

6270 

6275 

6280 

6285 

6290 

6295 

6300 

6305 

6310 

6315 

6320 

6525 

6330 

6335 

b 

6345 

6350 

6355 

6360 

6365 

6370 

6375 

6380 

6385 

6390 

6395 

6400 

6405 

6410 

6415 



PRINT "SR 1.5,2.PU Pi^ 40 , - 1 0 , DI ; LBq / ( MU/n . PR 0.5,1. SR 1,1.5; 
.PR .5,“1;LB)" 

ELSE 

PRINT "DI PA 38,-10.LB(T.PR .5,-l;LB£.PR .5,1;LB“T:PR -2.4,: 
PUiPR .5,-4;'* 

PRINT "LBuo.PR .5,1 ;LB)/K" 

END IF 

PRINT "CP 0,0 DI" 

Xlg*l.E+6 
Xug--1 . E+6 
Xal*50 
Yal*95 
Nrun*0 



BEEP 

INPUT “UANT TO PLOT DATA FROM A FILE (1 *Y , 0«N )? ” , Ok 
Xll-1 .E+6 



Xul*-1 .E+6 
Okp»0 

IF Ok*l THEN 
BEEP 

INPUT “ENTER THE NAME OF THE PLOT DATA F ILE " , D_f i le$ 
ASSIGN ©File TO D.file$ 

IF lconb<>0 THEN 6265 

Sx*0 

Sy*0 

Sx2«0 



Sxy=0 



Md«l 

BEEP 

INPUT "ENTER THE BEGINNING RUN NUMBER (DEF«l)",Md 

Npairs*9 

BEEP 

INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED (DEF»9)" 
Nrun*Nrun+Npa irs 



PRINTER IS 1 
BEEP 

PRINT USING "4X, 
PRINT USING "4X, 
PRINT USING "4X, 
PRINT USING "4X, 
PRINT USING "4X, 
PRINT USING "4X, 
INPUT Syn 



‘Select a symbol ^ “ 

" 1 Star 2 Plus sign""” 

3 Circle 4 Square""" 

" 5 Rombus""" 

" 6 Right-side-up triangle 

" 7 Up-side-down triangle""" 



bbtr 



Npairs 



INPUT "ENTER TUBE NUMBER FOR PLOTTING ( 0»TOP , l = SECOND , . . . ) " 

PRINTER IS 705 

IF Sym*l THEN PRINT "SM*" 

IF Sym=2 THEN PRINT "SM+" 

IF Sym«3 THEN PRINT "SMo" 

IF Md>l THEN 

FOR 1*1 TO (Md-1 ) 

ENTER 0Fi 1 e ; Xya( * ) 

NEXT I 
END IF 

FOR 1*1 TO Npairs 

ENTER @File;Xya(*) 

Ya*Xya( Itube*2 ) 

Xa*Xya ( 1 1 ube*2 + l ) 

Yc*LOG(Xa ) 



LB2iSF 
: PD PR 



, I tube 
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C w t 

6425 

£430 

6435 

5440 

6445 

6450 

6455 

6460 

6465 

6470 

6475 

6480 

6485 

6490 

6495 

6500 

6505 

6510 

6515 

6520 

6525 

6530 

6535 

6540 

6545 

6550 

6555 

6560 

6565 

6570 

6575 

6580 

6585 

6590 

6555 

6600 

6605 

6610 

6615 

6620 

6625 

6630 

6635 

6640 

6645 

6650 

6655 

6660 

6665 

6675 

6680 

6685 

6690 

6695 

6700 

6705 

6710 

6715 



xc=L06( Xa / } 

Sa=Sx+Xc 
5y=5y+Yc 
S^2=Sx2+>c' 2 
S;^y = S»y^Ac*Yc 
IF lht=0 THEN 
Xt*Ya 
Ya=va/Xe 
Xa = A t 

IF Xa/1 .E-^6>Xul THEN Xul=Xa/l.E + 6 
IF XaM.E + 6«'Xll THEN X 1 1 =Xa/ 1 . E-^6 
END IF 

IF lht=0 THEN 

X=(Xa*l .E-6-Xnin )*Sf ^ 

Y=( Ya* 1 . E-5-Ymin )*5fy 
END IF 

IF Iht=l THEN 

X = ( V a-Xn in ) * S f / 

Y*( YaM .E-6-Ynin )*Sfy 
END IF 

IF Iht=2 THEN 

X=' Xa/Ya-Xnin)*5f« 

Y=( Ya*l .E-3~Vnin )*Sfy 
END IF 

IF Y>100 OR Y<0 THEN 6585 
IF 5ym>3 THEN PRINT "Sn“ 

IF Syn'4 THEN PRINT "SR 1 . 4 , 2 . 4 *' 

PRINT ••PA",X,Y,"" 

IF Sym>3 THEN PRINT “SR 1. 2,1.6" 

IF Syn»4 THEN PRINT " UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , 8 , 4 , 0 ; “ 

IF Sym=5 THEN PRINT " UC3 , 0 , 99 , -3 , -6 , -3 , 6 , 3 , 6 , 3 . -6 ; " 
IF Syn=6 THEN PRINT " UC0 , 5 . 3 , 99 , 3 , “8 , “6 , 0 , 3 , 8 ; " 

IF Sym=7 THEN PRINT " UC0 , -5 . 3 , 99 , -3 , 8 , 6 , 0 , -3 , -8 ; " 
NEXT I 
BEEP 

INPUT "UANT TO LABEL (1 =Y , 0=N )'? ",I lb 1 
IF Ilbl=l THEN 

IF Syn>3 THEN PRINT -SM" 

IF Sym<4 THEN PRINT "SR 1.4, 2. 4" 

PRINT "PA" ,Xal ,Yal , 

IF Sym>3 THEN PRINT "SR 1.2, 1.6" 

IF Sym*=4 THEN PRINT "UC2 , 4 , 99 , 0 , *8 , -4 , 0 , 0 , 8 , 4 , 0 ; “ 

IF Syn*5 THEN PRINT " UC3 , 0 , 99 , -3 , -6 , -3 , 6 , 3 , 6 , 3 , -6 ; “ 
IF Syn*6 THEN PRINT “ UC0 , 5 . 3 , 99 , 3 , “8 , “6 , 0 , 3 , 8 ; " 

IF Syn*7 THEN PRINT " UC0 , -5 . 3 , 99 , -3 , 8 , 6 , 0 , -3 , -8 j " 
PRINT "SM" 

IF Sym<4 THEN PRINT "PR 2,0" 

PRINT “PR 2,-1.0;5R 1 . 0 , 1 . 8 ; LB " j D. f i 1 e$ ; " " 

Yal *Ya 1 -5 
BEEP 

IF Ias=l THEN 
BEEP 

INPUT "ENTER THE STRING ", Labe 1 $ 

PRINT "PR 2,0;SR 1 . 0 , 1 . 8 . LB " ; Labels ; " " 

GOTO 6665 
END IF 
END IF 
BEEP 

INPUT "WANT TO COMBINE ANOTHER FILE’ (1 =Y , 0-N ) " . I comb 
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6720 

6725 

6730 

6735 

6740 

6745 

6750 

6755 

6760 

6765 

6770 

6775 

6780 

6785 

6790 

6795 

6800 

6805 

6810 

6815 

6820 

6825 

6830 

6835 

6840 

6845 

6850 

6855 

6860 

6865 

6870 

6875 

6880 

6885 

6890 

6895 

6900 

6905 

6910 

6915 

6920 

6925 

6930 

6935 

6940 

6945 

6950 

6955 

6960 

6965 

6970 

6975 

6980 

6985 

6990 

6995 

7000 

7005 

7010 



ASSIGN ©File TO ♦ 

Xll«5 

Xul-45 

IF lconb<>0 THEN 6220 
BEEP 

INPUT "WANT TO PLOT A LEAST-SQUARES LINE ( 1 =Y , 0=N ) M 1 s 
IF Ils-l THEN 
BEEP 

INPUT "SELECT EXPONENT: 0*COMPUTE, 1 =0 . 75 ) " , lexp 
BEEP 

INPUT “SELECT CURUE TYPE ( 0 = SOLID, 1*DA5HED >M 1 1 

Ilt-Ilt+1 

PRINT "SM" 

IF Iexp-0 THEN 

Bb*( Nrun*Sxy-Sy*Sx )/ ( Nrun*Sx2-Sx ^^2 ) 

ELSE 

Bb*.75 
END IF 

Aa = ( SyBb^Sx )/Nrun 
Aa«EXP( Aa ) 

PRINTER IS 1 

PRINT USING "10X,""a = "\Z.4DE’‘;Aa 
PRINT USING "10X/'‘’n = 2 . 4DE " ; Bb 

PRINTER IS 705 
ln*0 

IF lht*0 THEN Xxstep=Xstep/40 
IF lht>0 THEN Xxstep=^Xstep/10 
FOR Xa-Xll TO Xul STEP Xxstep 
IF Xa>.99*Xnax THEN 6995 
IF Iht«l THEN Ya=Aa*Xa"Bb 

IF lht«0 THEN Ya«Aa''( 1/Bb )♦( Xa^l .E + 6 )"( ( Bb-1 >/Bb ) 
IF Iht-2 THEN Ya=Aa*Xa ( Bb-1 ) 

IF Iht-0 THEN 

Y»(Ya*l .E~3-Ynin)*Sfy 
X*( Xa-Xnin )*S f x 
END IF 

IF Iht-1 THEN 

Y«=(Ya^l.E-6-Ymin)^Sfy 
X*(Xa-Xnin)*Sfx 
END IF 

IF Iht-2 THEN 

Y*(YaM.E-3-Ynin>*Sfy 
X*( Xa-Xnin )*S f x 
END IF 

IF Y<0 THEN Y-0 
IF Y>100 THEN 6990 
IF Ilt«=l THEN 

PRINT “PA“,X,Y,'PD“ 

ELSE 

In*In+l 
Ir-In MOD lit 

IF Ir-1 THEN PRINT "PA‘’ , X . Y , "PD " 

IF Ir-0 THEN PRINT "PA“ , X , Y /’PU" 

END IF 
NEXT Xa 
PRINT “PU“ 
tNU Ih 
lconb«0 
GOTO 6185 
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7015 

7020 

7025 

7030 

7035 

7040 

7045 

7050 

7055 

70G0 

70G5 

7070 

7075 

7080 

7085 

7090 

7095 

7100 

7105 

7110 

7115 

7120 

7125 

7130 

7135 

7140 

7145 

7150 

7155 

71G0 

7165 

7170 

7175 

7180 

7185 

7190 

7195 

7200 

7205 

7210 

7215 

7220 

7225 

7230 

7235 

7240 

7245 

7250 

7255 

7260 

7265 

7270 

7275 

7280 

7285 

7290 

7295 

7300 

7305 



END IF 

PRINT "PU SM" 

BEEP 

INPUT "UANT TO PLOT NUSSELT LINE ( 1 = Y , 0=N ) 7 - , I np 

IF lnp=0 THEN 7125 

BEEP 

INPUT "ENTER TSAT (DEFAULT=18 DEG C)’\T5at 
Hf g*FNHfg{ Tsat ) 

Xil«5 

Xui-45 

FOR Xa^Xll TO Xul STEP Xstep/50 
Tf iln»Tsat-Xa* .5 
Kf-FNK(Tfiln) 

Rhof»FNRho(Tfiln ) 

Muf-FNMu(Tf iln) 

Ya*.728«( Kf "3*Rhof ''2*9.81*Hfg/(Muf*Do*Xa ) )“' .25 

X=(Xa-Xnin)*Sfx 

Y'(Ya*l .E-3-Ynin)*Sfy 

PRINT "PA" ,X, Y, "PD" 

NEXT Xa 

PRINT “PU PA 0,0" 

PRINT "PU PA 0,0 SP0" 

SUBEND 
SUB Plot3 

COM /Drl/ St ar , Syn , Icon 
COM /Fid/ Ift 
DIM C( 9 ) 

Fw= 1 

PRINTER IS 1 
BEEP 

PRINT USING "4X, ""Select Option X-Y Linits:""" 

PRINT USING "GX,""0 Use default values 

PRINT USING "6X,""1 Use new values 

INPUT Okd 
PRINTER IS 705 
IF Okd=0 THEN 
Xni n*0 
Yn 1 n=0 
Xmax= 15 
Yna^=15 
Xstep=3 
Ystep=3 
ELSE 
BEEP 

INPUT "ENTER MINIMUM AND MAXIMUM X-UALUES " , Xn i n , Xnax 
BEEP 

INPUT "ENTER MINIMUM AND MAXIMUM Y-UALUES " , Yni n , Ynax 
BEEP 

INPUT "ENTER STEP SIZE FOR X-AX IS “ , Xs t ep 
BEEP 

INPUT "ENTER STEP SIZE FOR Y-AXIS‘* , Ystep 
END IF 
BEEP 

PRINT "IN:SP1;IP 2300,1800,8300,6800:" 

PRINT “SC 0,100,0, 100;TL 2,0:" 

Sfx-100/(Xnax-Xmin) 

Sfy®100/(Ypnax-Ymin) 

BEEP 

lcg**0 
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7310 

7315 

7320 

7325 

7340 

7345 

7350 

7355 

7360 

7365 

7370 

7375 

7380 

7385 

7390 

7395 

7400 

7405 

7410 

7415 

7420 

7425 

7430 

7435 

7440 

7445 

7450 

7455 

7460 

7465 

7470 

7475 

7480 

7485 

7490 

7495 

7500 

7505 

7510 

7515 

7520 

7525 

7530 

7535 

7540 

7545 

7550 

7555 

7560 

7565 

7570 

7575 

7580 

7585 

7590 

7595 

7600 

7605 



INPUT "LIKE TO BY-PASS CAGE ( 1 =Y , 0 = N='DEFAULT )7 " , I eg 
IF Icg-1 THEN 7620 
PRINT "PU 0,0 PD" 

FOR Xa*^Xnin TO Xna^ STEP Xstep 
PRINT "PA";X,",0. XT. " 

NEXT Xa 

PRINT "PA 100,0. PU. " 

PRINT "PU PA 0,0 PD" 

POR Ya'Ynin TO Ymax STEP Ystep 
Y*( Ya-Ynin )*Sf y 
PRINT "PA 0, " . Y, "YT" 

NEXT Ya 

PRINT "PA 0,100 TL 0 2" 

FOR Xa*Xmin TO Xnax STEP Xstep 
X = ( Xa-Xni n )*Sf x 
PRINT "PA" iX, " ,100. XT" 

NEXT Xa 

PRINT "PA 100,100 PU PA 100,0 PD" 

FOR Ya=Ynin TO Ynax STEP Ystep 
Y = ^Ya-Yrniri)*Sfy 
PRINT "PD PA 100, " ,Y, "YT" 

NEXT Ya 

PRINT "FA 100, 100 PU" 

PRINT "PA 0,-2 SR 1.5.2" 

FOR Xa=Xnin TO Xnax STEP Astep 
X = ( Xa-Xr,in )*Sfx 
PRINT "PA",X,",0. " 

IF Xa<l AND XaO0 THEN PRINT "CP - 1 . 5 , - 1 ; LB0 ; PR - 1 , 0 ; LB " : Xa ; " 
IF Xa=0 THEN PRINT "CP -.5,-l;LB0" 

Xin=0 

IF Xa MOD 1*0 THEN Xin=l 

IF Xa> = 10 THEN PRINT "CP -2 , - 1 . LB " ; X a i " " 

IF Xa^l AND Xa<10 AND Xin*l THEN PRINT "CP - 1 . 25 , - 1 . LB “ ; Xa . " " 
IF Xa>l AND Xin=0 THEN PRINT "CF -2 , - 1 ; LB " ; Xa ; " " 

IF Xa*l THEN PRINT "CP -1,-1:LB1.0" 

NEXT Xa 

PRIN" "PU PA 0,0" 

Iht*2 ' modify 
FOR \a=Vpin TO Ynax STEP Ystep 
Y*( Ye-Yrun )*Sfy 
PRINT "PA 0, " i Y , " " 

IF lht*0 AND Ya.'0 THEN PRINT "PR 2,0" 

IF Ya<l AND Ya<>0 THEN PRINT "CP -4 , - . 25 . LB0 . PR -2 , 0 . LB " i Ya ; " 
IF Ya*0 THEN PRINT "CF -2,-.25.LB0" 

IF Ya>l AND Iht<2 THEN PRINT "CP -5 , - . 25 ; LB " . Ya . " " 

IF Ya>9 AND Iht=2 THEN PRINT "CP -4 , - . 25 ; LB " ; Ya ; " " 

IF Ya>0 AND Ya<10 AND Iht = 2 THEN PRINT "CP -3 , - . 25 ; LB “ ; Ya ; " 
IF Ya*l THEN PRINT "CP -4 , - . 25 ; LB 1 . 0 " 

NEXT Ya 
IF OLd*0 THEN 
XlabelS*"X" 

Ylabei$="Y" 

ELSE 

BEEP 

INPUT "ENTER X-LABEL " , X 1 abe 1 $ 

BEEP 

INPUT "ENTER Y-LABEL " , Y1 abel $ 

END IF 

PRINT "SR 1.5,2 PU PA 50,-10 CP" ; -LEN( Xlabel$ )/2 ; "0 ; LB ‘ ; XlabelS; " 
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7610 

7615 

7620 

7625 

7630 

7635 

7640 

7645 

7650 

7655 

/ClDW 

7665 

7670 

7675 

7680 

7685 

7690 

7695 

7700 

7705 

7710 

7715 

7720 

7725 

7730 

7735 

7740 

7745 

7750 

7755 

7760 

7765 

7770 

7775 

7780 

7785 

7790 

7795 

7600 

7805 

7810 

7615 

7820 

7825 

7630 

7835 

7840 

7845 

7850 

7855 

7860 

7865 

7870 

7875 

7880 

7885 

7890 

7895 

7900 



PRINT "PA -11,50 CP 0," .-LEN(Ylabel$)/2*5/6; "DI 0 , 1 . LB " ; Ylabel $ j “ 

PRINT "CP 0,0 or* 

Nrun«»0 

BEEP 

INPUT “WANT TO PLOT DATA FROM A FILE (1 -Y , 0*N >7 " , Ok 
Okp = 0 

IF Ok*l THEN 
BEEP 

INPUT "ENTER THE NAME OF THE PLOT DATA F I LE " , D_f 1 1 e$ 

ASSIGN ©File TO D_file$ 

It- ItitN 

Sx = 0 
Sy*0 
Sx2 = 0 
Sxy = 0 
Md*l 



BEEF 

INPUT "ENTER THE BEGINNING RUN NUMBER (DEF = n",Md 

Npa 1 rs=9 

BEEP 

INPUT "ENTER THE NUMBER OF X-^ PAIRS STORED ( DEF = 9 ) " , Npa 1 r 5 
BEEP 

INPUT "SELECT TUBE NUMBER ( 0=TOP , 1 =SEC0ND )",Itube 

Nrun=Nrun+Npairs 



PRINTER IS 1 
BEEP 

PRINT USING "4X,""Sei 
PRINT USING "4X,"" 1 
PRINT USING "4X , " " 3 
PRINT USING "4X,“" 5 
PRINT USING "4X,"" 6 
PRINT USING "4X,"“ 7 
INPUT Syn 
PRINTER IS 705 



ec t a synbo 1 ; " " " 

Star 2 Plus sign 

Circle 4 Square 

Rombus " " " 

R 1 ght -5 1 de-up triangle"" 
Up-s 1 de-down triangle 



IF Syn=l THEN PRINT "SM*" 
IF Syn*2 THEN PRINT "SM+" 
IF Syn=^3 THEN PRINT “SMo" 
Md-Itube*9 



IF Md/1 THEN 



FOR I»0 TO ( Md-1 > 
ENTER @File.X,Y 
NEXT I 
END IF 

FOP 1=1 TO Npairs 

ENTER @File,Xa,Ya 

Sx^Sx'^'Xa 

Sy*Sy+Ya 

Sx2*Sx2+Xa"2 

Sxy=Sxy+Xa*Ya 

X = ( Xa-Xrun )*Sf a 



Y = ( Ya-Ynin )*Sf y 
IF Y>100 OR Y<0 THEN 7910 
IF Syrn>3 THEN PRINT 
IF Syrn<4 THEN PRINT 
PRINT "PA",X,Y,"" 

IF Syrn>3 THEN PRINT 
IF Syrn*4 THEN PRINT 
IF Syn*5 THEN PRINT 
IF Syn-6 THEN PRINT 



SM" 

SR 1. 4,2.4 " 

SR 1.2, 1.6" 

UC2,4,99,0,-8,-4,0,0,8,4,0j " 
UC3,0,99,-3,-6,-3,G,3,G,3,-6; " 
UC0,5.3,99,3,-8,-6,0,3,8; " 
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7905 

7910 

7915 

79Z0 

79Z5 

7930 

7935 

7940 

7945 

7950 

7955 

7960 

7965 

7970 

7975 

7980 

7985 

7995 

8000 

8005 

8010 

8015 

8020 

8025 

8030 

8035 

8040 

8045 

8050 

8055 

8060 

8065 

8070 

8075 

8080 

8085 

8090 

8095 

8100 

8105 

8110 

6115 

8120 

6125 

8130 

8135 

8140 

6145 

8150 

8155 

8160 

8165 

8170 

8175 

8180 

8185 

8190 

8195 

8200 

8205 



IF 5yn«=7 THEN PRINT “ UC0 , -5 . 3 , 99 , -3 , 8 , 6 , 0 , “3 , -6 ; " 

NEXT I 
BEEF 

INPUT ’’WANT TO LABEL < 1 = Y , 0 = N )7 ",1 1 b 1 
IF Ilbl=l THEN 

IF 5yn>3 THEN PRINT •'SM" 

IF 5yn^4 THEN PRINT "SR 1.4, 2. 4" 

PRINT "PA" ,Xal ,Yal , " “ 

IF 5yn>3 THEN PRINT "SR 1.2, 1.6" 

IF Syn*4 THEN PRINT “ UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , 8 , 4 , 0 ; " 

IF Sym-5 THEN PRINT " UC3 , 0 , 9 9 , -3 , -6 , -3 , 6 , 3 , 6 , 3 , -6 . ‘ 

IF Syn=6 THEN PRINT “000,5.3,99,3,-8,-6,0,3,8;" 

IF Syn-7 THEN PRINT ” UC0 , -5 . 3 , 99 , -3 , 8 , 6 , 0 , -3 , -8 ; " 

PRINT “SM" 

IF 5yn<4 THEN PRINT "PR 2,0" 

PRINT "PR 2,-1.0;SR 1 . 0 , 1 . 8 ; LB " ; D_f i 1 e$ : " " 

Yal =Yal -5 

INPUT "UANT TO ADD ANOTHER STRING ( 1 =Y , 0-N , las 

IF las^l THEN 
BEEP 

INPUT "ENTER THE 5TR I NG " , Labe 1 $ 

PRINT "FR 2,0:SR 1 . 0 , 1 . 8 . LB " . Labe 1 $ ; " " 

GOTO 7990 
END IF 
END IF 
BEEP 

INPUT "WANT TO COMBINE ANOTHER FIlE'? ( 1 =Y , 0 = N ) " , Iconb 
ASSIGN @File TO * 

IF IconbOO THEN 7645 

Ils = l 

BEEP 

INPUT "UANT TO PLOT A LEAST-SQUARES LINE ( 1 =DEF = YES , 0=NO ) " , 1 1 s 
IF Il5=l THEN 
BEEP 

INPUT "SELECT CURUE TYPE < 0-SOL I D , 1 -DASHED 1 1 1 

Ilt-Ilt+1 

PRINT "SM" 

IF lexp=0 THEN 

Bb=^Nrun*S>:y-Sy*S-* )/(Nrun*SA2-S>, '2 ) 

ELSE 

Bb=.”5 
END IF 

Aa-< Sy-Bb*Sx )/Nrun 
PRINTER IS 1 

PRINT USING "10X,""a = " " , MZ . 3DE " : Aa 
PRINT USING "10X.""n - " " , MZ . 3DE " ; Bb 
PRINTER IS 705 
In-0 

FOR Xa=Xnin TO Xnax STEP (Xnax-Xnin) 

Ya=Aa + Xa *Bb 

Y=< Ya-Yrun )*Sf y 

X = < Xa-Xmin )*Sfx 

IF Y<0 THEN Y=0 

IF Y>100 THEN GOTO 8220 

IF Ilt=l THEN 

PRINT "PA" ,X, Y, "PD" 

ELSE 

I n=In+ 1 
Ir-In MOD lit 

IF Ir=l THEN PRINT " PA " , X , Y , " PD " 
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6Z10 

6:!5 

622C 

8225 

6230 

8235 

8240 

6245 

8250 

8255 

8260 



IF lr*0 THEN PRINT "PA" , X , Y , “PU" 
END IF 
NEXT Xa 
PRINT “PU" 

END IF 
I conb=0 
GOTO 7820 
END IF 

PRINT "PU PA 0,0" 

PRINT "PU PA 0,0 SP0" 

SUBEND 
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